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1.  Introduction 


The  study  of  particle  deposition  from  aerosol  flow  onto  surfaces  in  enclosed  spaces  has  attracted  con¬ 
siderable  attention  in  the  past  few  decades  owing  to  its  significance  in  a  great  variety  of  technological 
applications,  such  as  microcontamination  control  in  the  semiconductor  industry,1,  aerosol  production 
in  chemical  reactors/  indoor  air  pollution,44  and  soiling  of  artworks,7  as  well  as  nuclear  reactor 
safety  analysis.8  9  Both  desirable  and  undesirable  outcomes  may  result  from  the  processes  of  particle 
deposition.  For  instance,  the  thin-film  coating  techniques  used  to  generate  unique  surface  properties 
are  an  application  of  the  desirable  outcome.10 11  In  contrast,  there  are  many  examples  of  undesirable 
particle  deposition,  which  causes  surface  contamination  or  even  material  damage.  One  of  the  well- 
known  examples  is  the  reduced  yield  of  integrated  circuits  in  the  wafer  manufacturing  process.  Par¬ 
ticulate  contamination  on  critical  surfaces  of  a  space  telescope  leading  to  optical  obscuration  due  to 
light  scattering  and  hence  degrading  performance  is  another  unwanted  consequence.1 

To  develop  effective  strategies  to  control  either  desirable  or  undesirable  deposition  of  particles,  a 
sound  understanding  of  the  underlying  particle  transport  mechanisms  and  the  associated  physical 
factors  influencing  the  processes  is  required  to  gain  insights  into  particle  deposition  behaviors.  A 
significant  body  of  scientific  research  has  been  devoted  to  studying  the  dynamic  behavior  of  particles 
under  well-defined  configurations  and  air  flow  conditions — examples  include  aerosol-laden  air  flow¬ 
ing  through  fibrous  filters, n  14  particle  deposition  from  turbulent  channel  flow,15’16  and  particle  depo¬ 
sition  from  laminar,  vertical  flow  onto  a  circular  plate.2 17  19  In  contrast,  not  many  published  studies 
consider  particle  deposition  onto  surfaces  in  a  three-dimensional  confined  space,  probably  due,  in 
part,  to  the  complexity  of  airflow  patterns  associated  with  enclosures  of  various  geometrical  shapes. 

The  goal  of  this  report  is  to  provide  a  state-of-the-art  literature  review  concerning  particle  deposition 
onto  surfaces  within  an  enclosed  volume.  The  knowledge  gleaned  from  the  available  scientific  publi¬ 
cations  presented  here  can  be  applied  to  an  enclosure  scale  as  small  as  a  spacecraft  payload  cavity,  or 
as  large  as  a  clean  room  in  a  manufacturing/processing  facility.  The  surfaces  of  interest  for  particles 
to  deposit  on  could  be  walls,  the  floor,  the  ceiling,  or  any  contamination-sensitive  surfaces  within  an 
enclosed  volume  relevant  to  various  engineering  applications. 

The  report  begins  with  the  characterization  of  airflow  adjacent  to  a  surface  and  how  it  relates  to  the 
core  region  of  an  enclosure.  Particle  deposition  flux,  a  parameter  used  to  evaluate  the  rate  of  deposi¬ 
tion  onto  surfaces,  is  briefly  introduced  as  the  concept  will  be  mentioned  throughout  the  chapter. 

Next,  the  fundamental  physical  processes  governing  the  transport  of  particles  are  discussed,  and  the 
definition  of  the  parameters  commonly  reported  in  the  literature  for  characterizing  deposition  are 
introduced,  along  with  the  available  experimental  techniques  for  their  measurement.  The  dependence 
of  particle  deposition  rate  on  particle  size,  airflow,  and  surface  characteristics  are  examined  from  the 
existing  published  experimental  investigations  relevant  to  enclosures.  Lastly,  modeling  developments 
for  predicting  the  rate  of  particle  deposition  and  associated  experimental  validations  are  presented. 
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2.  Background 


Knowledge  of  particle  transport  and  deposition  within  an  enclosed  volume  has  been  applied  in 
numerous  practical  scenarios.  In  the  beginning  of  this  section,  special  attention  is  given  to  the  rele¬ 
vance  of  particle  deposition  to  contamination  control  in  the  context  of  space  applications. 

Since  aerosol  particles  are  transported  to  the  vicinity  of  surfaces  by  air  currents,  understanding  the 
nature  of  air  flow  and  its  near-surface  characteristics  is  the  first  step  to  gain  insights  into  the  processes 
of  particle  deposition  onto  surfaces.  In  addition,  a  quantitative  approach  using  the  concept  of  flux  and 
the  underlying  physical  mechanisms  that  govern  particle  deposition  are  briefly  described. 


2.1  Relevance  to  Contamination  Control  in  Space  Applications 

In  space  applications,  the  concern  of  particle  deposition  on  contamination-sensitive  surfaces  arises 
from  the  undesirable  effects  of  surface  obscuration.  For  one  thing,  the  presence  of  particles  on  optical 
reflective  surfaces,  e.g.,  mirrors  and  focal  planes,  interferes  with  the  proper  transmission  of  light  to 
the  next  optical  components,  hence  reducing  the  signal  strength.  Particles  residing  on  an  absorptive 
surface,  such  as  baffles  inside  a  telescope  assembly,  causes  light  scattering  , which,  in  turn,  over¬ 
whelms  the  signal  and  ultimately  may  even  damage  other  optical  components.  The  presence  particles 
on  thermal  control  surfaces  causes  alterations  of  solar  absorbance  and/or  emissivity,  which  leads  to  an 
altered  equilibrium  temperature,  and  may  further  deteriorate  the  thermal  control  function. 

In  light  of  the  adverse  effects  due  to  particulate  contamination,  the  aim  of  contamination  control  is  to 
prevent  performance  degradation  by  minimizing  the  deposition  of  particles  on  critical  and  sensitive 
surfaces  of  spacecraft  components.  Meeting  this  objective  is  particularly  important  for  high-profile 
remote  sensing  spacecraft,  such  as  the  Hubble  Space  Telescope.  As  a  consequence,  tremendous 
efforts  and  resources  are  dedicated  to  contamination  control  during  various  phases  of  spacecraft  proc¬ 
essing:  design,  manufacturing,  assembly,  testing,  integration,  storage,  shipping,  launch  site  ground 
activities,  and  in-flight  operations. 

The  development  of  effective  mitigation  strategies  to  minimize  particle  contamination  requires 
knowledge  of  particle  transport  and  deposition  as  well  as  the  associated  physical  factors  affecting  the 
processes.  Spacecraft  surface  contamination  as  a  result  of  particle  fallout  may  take  place  at  any  level 
of  engineering  practice,  for  instance,  component  fabrication  in  a  cleanroom ,  system  assembly  and 
testing  within  a  payload  cavity ,  and  spacecraft  integration  inside  payload  fairing.  Note  that  the  terms 
in  italics  refer  to  the  corresponding  enclosed  volumes  under  the  example  scenario  considered;  thus, 
one  can  see  that  the  study  of  particle  deposition  onto  surfaces  in  enclosures  is  of  strong  relevance  to  a 
variety  of  circumstances  in  space  applications.  The  insights  from  studying  the  physical  processes  that 
influence  the  rate  of  particle  deposition  within  an  enclosure  are  expected  to  provide  a  strong  scientific 
foundation  to  benefit  various  aspects  of  aerospace  contamination  control  needs. 
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2.2  Nature  of  Air  Flow  and  its  Near-Surface  Characteristics 

Inside  an  enclosed  volume,  air  is  mixed  through  two  mechanisms:  natural  convection  and  turbulent 
mixing  (forced  convection).  Natural  convection  is  driven  by  temperature  gradients,  i.e.,  temperature 
difference  (AT)  across  the  space  inside  the  enclosure,  induced  by  heat  transfer  at  surfaces.  Turbulent 
mixing  or  forced  convection,  on  the  other  hand,  is  associated  with  external  energy  input  into  a  system, 
such  as  fan  mixing  or  fluid  flushing.  In  an  enclosure  without  the  mechanism  of  forced  convection, 
natural  convection  becomes  the  only  cause  for  fluid  mixing.  Small  temperature  variations  within  the 
enclosure  may  induce  convective  currents,  which,  in  turn,  cause  aerosol  concentration  to  become 
mixed  throughout  the  volume.  It  is  estimated  that  the  convective  flow  velocity  can  reach  approxi¬ 
mately  1  cm/s  in  a  1-m-high  chamber  when  the  floor  or  wall  surfaces  are  warmer  than  the  air  by 
0.01  °C.20  Study  also  has  shown  that  a  AT  of  0.1  °C  can  induce  homogeneous  mixing  for  particles  up 
to  15  |im  in  diameter.21  On  the  other  hand,  buoyancy  can  also  impede  mixing.  For  instance,  a  warm 
ceiling  and  cool  floor  will  promote  stratification  of  air  masses  within  an  enclosure.  Some  empirical 
evidence  about  near-surface  airflows  in  rooms  can  be  found  in  reference  [22] 

One  widely  accepted  approach  to  model  air  flow  in  enclosures  is  to  assume  that  the  air  in  the  core 
region  is  homogeneously  and  isotropically  turbulent,  behaving  like  an  ideal  nonviscous  fluid.  Adjacent 
to  the  interior  surface,  the  air  is  assumed  to  behave  as  a  viscous  fluid  within  a  thin  layer,  which  is  also 
known  as  the  boundary  layer.  Inside  the  momentum  boundary  layer,  the  fluid  velocity  drops  sharply  to 
zero  at  the  surface  from  its  mainstream  peak  value.  The  thickness  of  the  momentum  boundary  layer 
depends  on  the  momentum  diffusivity  of  the  fluid.  For  instance,  the  typical  thickness  is  about  1  to  2  cm 
for  room  air  motion,6  given  air  kinematic  viscosity  of  0. 1 5  cm2  s-1  at  293K  and  1  atm. 

The  boundary  layer  concept  can  also  be  used  to  capture  the  near-surface  characteristics  in  terms  of 
thermal  and  contaminant  concentration  profiles.  For  air,  the  thickness  of  the  thermal  boundary  layer 
is  comparable  to  that  of  the  momentum  boundary  layer  because  the  thermal  diffusivity  (-0.2  cm2  s~! 
for  air)  is  of  the  same  order  of  magnitude  as  the  momentum  diffusivity.  For  gaseous  contaminants 
with  molecular  weights  close  to  those  of  N2  and  02,  the  thickness  of  the  concentration  boundary  layer 
would  be  similar  to  that  of  the  momentum  and  thermal  boundary  layer  because  the  molecular  diffu- 
sivities  are  comparable  to  momentum  and  thermal  diffusivities.  Airborne  particles,  however,  have 
small  diffusivities  (or  diffusion  coefficients)  compared  to  molecular  contaminants,  which  leads  to 
much  thinner*  particle  concentration  boundary  layers  than  the  momentum  and  thermal  boundary  lay¬ 
ers.  Okuyama  et  al.2^  estimated,  from  their  experiments  in  a  2.6-L  non-stirred  cylindrical  tank,  the 
particle  concentration  boundary  layer  thickness  to  be  approximately  0.3  cm  when  the  particle  diffu¬ 
sivity  is  \0~'  cm2  s“\  Considerable  information  about  particle  concentration  boundary  layer  thickness 
can  be  found  in  references  [24]  and  [25]. 

In  addition  to  aerosol  diffusivities,  the  nature  and  the  intensity  of  the  near-surface  air  flow  also  play  a 
role  in  determining  the  thickness  of  the  concentration  boundary  layer.  Scale  analysis26  that  is  used  to 
approximate  particle  boundary  layer  thickness  adjacent  to  surfaces  indicates  that,  given  the  same  dif¬ 
fusivity,  the  boundary  layers  are  thinner  when  the  fluid  outside  them  is  fast  moving,  and  thicker  when 
the  fluid  moves  slowly. 


As  will  be  shown  later  in  Subsection  6.3.2,  the  implication  of  a  thinner  concentration  boundary  layer  is  that  small-scale 
surface  roughness  can  play  a  significant  role  in  affecting  particle  transport  across  the  boundary  layer,  but  have  negligible 
effects  on  momentum  and  heat  transfer. 
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In  general,  the  air  within  an  enclosure  of  arbitrary  shape  can  be  visualized  to  consist  of  two  parts:  a 
core  zone  where  the  air  is  well  mixed,  and  a  thin  quiescent  boundary  layer  adjacent  to  the  inside  sur¬ 
face  of  the  enclosure  where  little  air  motion  exists  in  the  direction  perpendicular  to  the  surface  (see 
Figure  1).  In  the  core  zone  the  particle  concentration  is  often  assumed  to  be  spatially  uniform  due  to 
well-mixed  air  flows,  and  large-scale  turbulent  (or  eddy)  diffusion  is  responsible  for  bringing  aerosols 
to  the  vicinity  of  the  surface  for  subsequent  deposition/  Upon  arrival  at  the  boundary  layer,  aerosols 
may  migrate  through  the  thin  layer  to  the  surface  by  mechanisms  such  as  Brownian/turbulent  diffu¬ 
sion,  gravitational  settling,  thermophoresis,  inertial  drift,  and  electrostatic  attraction.  These  transport 
processes,  as  will  be  discussed  in  Section  3,  commonly  control  the  rate  of  particle  deposition  onto 
surfaces  in  enclosures. 

In  summary,  the  boundary  layer  concept  used  to  characterize  aerosol  transport  processes  from  main¬ 
stream  fluid  onto  surfaces  has  been  applied  with  great  success  in  a  large  variety  of  air  flow  scenarios, 
such  as  external  flows  around  cylinders  or  spheres,  laminar  and  turbulent  pipe  flows,  and  enclosures 
with  quiescent  or  turbulent  flows.  Theoretical  developments  based  on  this  boundary  layer  concept  to 
evaluate  the  extent  of  particle  deposition  within  an  enclosed  volume  along  with  the  associated 
experimental  validations  will  be  discussed  in  more  details  in  Subsection  7.1. 


2.3  Evaluating  Particle  Deposition:  Deposition  Flux 

The  rate  of  particle  transport  onto  a  surface  is  commonly  quantified  in  terms  of  a  deposition  flux ,  with 
the  dimensions  of  mass  (or  number)  per  unit  time  per  unit  area.  In  the  event  that  transport  through  the 
boundary  layer  is  dominated  by  the  combined  effects  of  turbulent  and  Brownian  diffusion,  the  flux  of 
particles  onto  a  smooth,  isothermal,  and  electrically  neutral  vertical  surface  can  be  characterized  by 
the  following  one-dimensional  steady-state  continuity  equation. 

Flux  =  -(e_  +  ,  (1) 

p  dy 


1 
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quiescent 
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boundary  layer 


Figure  1.  Schematic  of  the  core  zone  and  the  boundary  layer  in  enclosures  of  different  geometrical  shapes. 


Also  known  as  "homogeneous  turbulence  model,”  and  its  theoretical  representations  will  be  addressed  in  Subsection  7.1. 
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where  £p  and  D  are  the  turbulent  (or  eddy)  and  Brownian  diffusivities  for  aerosols,  respectively;  C  is 
the  aerosol  concentration  in  air;  and  y  is  the  distance  from  the  surface. 

Equation  (1)  can  be  viewed  as  a  modified  form  of  Fick’s  first  law  of  diffusion ,  which  describes  the 
linear  relationship  between  the  flux  of  aerosols  and  the  concentration  gradient  dC/dy,  with  the  term 
(£p+  D)  being  the  proportionality  constant.  Physically,  this  suggests  that  the  net  transport  of  particles, 
due  to  diffusion  only,  always  takes  place  from  regions  of  high  to  low  particle  concentration. 

Equation  (1)  has  long  been  used  to  evaluate  particle  deposition  flux.  Since  the  air  flow  outside  the 
boundary  layer  is  assumed  to  be  homogeneously  mixed,  a  particle  concentration  gradient  only  exists 
within  the  boundary  layer.  Other  assumptions  for  this  theoretical  representation  include:  (1)  particles 
are  completely  retained  once  they  collide  onto  the  surface,  i.e.,  the  surface  acts  as  a  perfect  sink  for 
particle  deposition;  and  (2)  no  mechanisms  of  coagulation,  condensation,  and  evaporation  are 
involved  during  the  transport  process  (i.e.,  there  are  no  sources  or  sinks  for  particles  within  the 
boundary  layer);  thus,  the  particle  flux  is  constant  throughout  the  particle  concentration  boundary 
layer. 

In  the  presence  of  particle  transport  mechanisms  other  than  diffusion,  the  particle  flux  can  be  evalu¬ 
ated  by  adding  terms  to  account  for  external  force  fields  acting  on  the  particles,  and  that  gives 

Flux  =  -(e  +  D)^-  -  veC  ,  (2) 

dy 


where  ve  is  the  particle  steady-state  drift  velocity  under  the  action  of  external  body  forces  counter¬ 
acted  by  fluid  drag.  The  external  force  fields  of  most  interest  in  particle  deposition  processes  include 
gravitational,  electrical,  and  thermal  (generated  by  a  temperature  gradient  in  the  fluid),  as  will  be  pre¬ 
sented  in  the  next  section. 
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3.  Mechanisms  of  Particle  Transport 


Airborne  particles  are  transported  onto  surfaces  through  a  variety  of  physical  processes,  so-called 
deposition  mechanisms.  The  fundamental  physics  behind  particle  transport  or  movement  from  one 
point  to  another  is  universal,  regardless  the  configurations  of  the  system,  for  instance,  dust  accumula¬ 
tion  in  ventilation  ducts,  gas  cleaning  by  particle  collection  on  fibrous  filters,  and  scavenging  of  par¬ 
ticulate  matter  from  the  atmosphere.  Below,  the  underlying  physical  processes  governing  particle 
motion  in  air  are  addressed. 


3.1  Brownian  Diffusion 

Brownian  diffusion  is  the  characteristic  random  wiggling  motion  of  small  airborne  particles  in  still 
air,  resulting  from  constant  bombardment  by  surrounding  gas  molecules.  Such  irregular  motions  of 
pollen  grains  in  water  were  first  observed  by  the  botanist  Robert  Brown  in  1827,  and  later  similar 
phenomena  were  found  for  small  smoke  particles  in  air.  In  the  early  twentieth  century,  the  relation¬ 
ships  characterizing  Brownian  diffusion  based  on  kinetic  theory  of  gases  were  first  derived  by  Ein¬ 
stein  in  the  early  1900s,  and  later  verified  through  experiments.20 

The  Brownian  diffusivity,  D  (or  diffusion  coefficient),  which  relates  the  gas  properties  and  the  parti¬ 
cles  through  fluid  drag,  can  be  evaluated  by  the  Stokes-Einstein  expression.  7  K 


kTCc 
37tpdp  ’ 


(3) 


where  k  is  Boltzmann’s  constant  (1 .38  x  10  '23  J  K"1),  T  is  the  absolute  temperature  in  K,  1 1  is  dynamic 
viscosity  of  air,  dp  is  particle  diameter,  and  Cc  is  a  slip  correction  factor  for  small  particles  (see  Sub¬ 
section  3.3  for  details).  The  value  of  D  depends  on  the  particle  size  and  fluid  properties.  For 
instance,  the  diffusivity  of  a  0.01-pm  aerosol  particle  is  20,000  times  higher  than  that  of  a  10  pm 
aerosol  particle.  The  larger  the  value  of  D,  the  more  rapid  the  mass  transfer  process  to  drive  particles 
moving  from  regions  of  high  to  low  concentration.  Brownian  diffusion  is  the  dominant  particle  depo¬ 
sition  mechanism  for  small  particles  (<0.1  pm)  over  short  distances. 


3.2  Turbulent  Diffusion 

For  particle  transport  on  greater  physical  scales  (i.e.,  larger  distance),  turbulent  diffusion  is  far  more 
effective  than  the  thermal-based  Brownian  diffusion.  Mass  transfer  of  particles  from  one  point  to 
another  by  turbulent  (or  eddy/convective)  diffusion  is  analogous  to  momentum  transfer  by  fluctuating 
velocity  components  in  turbulent  flows.  Unlike  momentum  transfer  inside  the  thin  stagnant  layer  of 
air  adjacent  to  a  surface  where  viscous  forces  are  dominant  and  turbulence  fluctuations  are  negligible, 
mass  transfer  of  particles  is  attributed  to  the  combination  of  turbulent  and  Brownian  diffusion.  A 
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general  expression  for  the  particle  flux  to  the  surface  due  to  both  turbulent  and  Brownian  diffusion 
was  described  in  Eq.  (1),  with  the  derivations  shown  in  reference  [27]. 

Intuitively,  the  particle  turbulent  diffusivity,  £p,  inside  the  boundary  layer  is  expected  to  vary  with 
distance  from  the  surface,  y,  owing  to  the  physical  constraints  imposed  by  the  surface.  In  fact,  it  is 
extremely  difficult  to  determine  the  value  of  £p  explicitly  because  it  is  strongly  associated  with  the 
degree  of  turbulence  with  respect  to  the  air  flow  structure  in  the  system,  as  well  as  with  the  size  of 
particles  (e.g.,  heavy  particles  cannot  faithfully  follow  the  fluid  eddies).  Hinze29  has  shown  that  over 
a  sufficiently  long  time  scale,  the  particle  turbulent  diffusivity  equals  the  fluid  turbulent  viscosity  by 
solving  the  equation  of  motion  for  a  particle  in  a  homogeneous  turbulent  flow  field.  A  similar  argu¬ 
ment  made  by  Fuchs20  is  that  the  motion  of  larger  particles  is  more  persistent  due  to  their  larger  mass, 
which,  in  turn,  contributes  to  the  nearly  identical  average  distance  traveled  by  large  and  small  parti¬ 
cles  over  a  long  time  limit.  Additionally,  numerical  simulation  of  turbulent  diffusion  of  particles  also 
indicated  that  particle  and  fluid  turbulent  diffusion  are  comparable  for  particles  smaller  than  approxi¬ 
mately  170  |im.™  Therefore,  it  is  a  fair  statement  that  the  turbulent  diffusivity  for  particles  can  be 
reasonably  approximated  as  the  same  order  of  magnitude  as  that  of  gas  fluid  (-KT1  cm2  s_I),  which  is 
far  greater  than  the  particle  Brownian  diffusivities.  In  other  words,  particles  tend  to  migrate  closer  to 
the  surface  by  turbulent  diffusion  before  Brownian  diffusion  becomes  important. 

In  essence,  particle  turbulent  diffusivity  gradually  diminishes  in  the  particle  concentration  boundary 
layer  to  zero  at  the  surface.  Throughout  the  core  of  an  enclosure,  turbulent  diffusivities  can  vary  over 
large  ranges  in  accordance  with  the  intensity  of  air  motion. 


3.3  Drag  Force 

Fluid  drag  on  a  particle  is  the  resistance  force  exerted  by  the  surrounding  fluid  when  there  is  relative 
motion  between  the  particle  and  the  fluid.  In  other  words,  the  drag  force  is  always  present  as  long  as 
the  particle  is  not  travelling  in  a  vacuum.  The  net  effect  of  the  drag  force  is  to  reduce  the  acceleration 
of  particles. 

The  drag  force,  FD,  on  a  spherical  particle  can  be  characterized  by  Stokes’s  law: 

Fd  =  3njiVdp, 


where  V  is  the  particle  speed  relative  to  the  local  fluid  speed.  Stokes’s  law  is  derived  from  the  solu¬ 
tion  of  the  Navier-Stokes  equations,  assuming*  that  inertial  forces  are  negligible  compared  to  viscous 
forces  and  that  the  fluid  speed  at  the  particle’s  surface  is  the  same  as  the  particle’s.  Owing  to  the  low 
velocities  and  small  particle  sizes  involved,  most  aerosol  motions  occur  at  low  particle  Reynolds 
numbers  (Rep  «  1),  where  Rep  =  dpV/v,  expressing  the  ratio  of  inertial  to  viscous  forces  on  particles. 


Other  assumptions  include:  incompressible  fluid,  no  walls  or  other  particles  nearby,  constant  motion,  zero  fluid  velocity 
at  the  particle’s  surface,  and  rigid  spherical  particles.  These  assumptions  work  well  in  most  cases  for  aerosols.  See 
Hinds28  for  more  detailed  discussions. 


(4) 
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Table  1  shows  the  calculated  Rep  for  spherical  particles  of  different  sizes  falling  at  their  steady-state 
settling  velocity  in  air  due  to  gravity.  By  looking  at  Rep,  Stokes’s  law  clearly  holds  true  for  particles 
smaller  than  20  pm.  For  aerosols  that  are  relatively  large  and  move  through  a  fluid  rapidly,  the  iner¬ 
tial  forces  become  dominant  compared  to  viscous  forces;  thus,  the  drag  force  exerted  to  the  particle  is 
calculated  as 


FD  = 


f 

V 


2 


where  Cd  is  the  drag  coefficient,  given  by  reference  [31] 


(5) 


r 


cd  =  < 


24 

Re 

O  Q 

1  +  — Re+  — Re2  ln(2Re) 
16  160 

—  (l+0.15Rea687) 

Re v  ’ 

0.44 


24 

Re 


Re  «  1  (Stokes’s  law) 

0.1  <  Re  <  2 

2  <  Re  <  500 
500  <  Re  <  2xl05 


(6) 


When  the  size  of  an  aerosol  particle  approaches  the  mean  free  path  of  gas  molecules  (0.066  |im  for  air 
at  1  atm  and  20°C),  the  assumption  of  zero  relative  fluid  velocity  right  at  the  particle  surface  fails, 
which  could  lead  to  significant  errors.  In  this  case,  the  actual  drag  force  is  smaller  than  predicted  by 
Stokes’s  law.  A  correction  parameter,  the  Cunningham  correction  factor,  Cc,  is  introduced  to  correct 
for  this  “slip”  phenomenon  for  small  particles,  and  the  corrected  Stokes’s  law  becomes 


Fd  = 


37ipVdp 

Cc 


(7) 


Table  1 .  Particle  Reynolds  Numbers  (Rep)  Calculated  for  Parti¬ 
cles  of  Different  Diameters  Falling  at  Their  Terminal 
Settling  Velocities  in  Air  at  20°C,  Pressure  =  1  atm,  and 
Gravitational  Acceleration  =  980  cm  s*2 


Particle  diameter,  pm 

Rep 

0.1 

5.8x10  9 

1 

2.3x1  O’6 

10 

2.0x1  O'3 

20 

1.6x10  2 

50 

0.25 

100 

2 

300 

20 

9 


(8) 


Cc 


(  dDVI 

2.514  +  0.800exp 

-0.55-2- 

where  X  is  the  mean  free  path  of  gas  molecules.  The  empirical  formula  of  Eq.  (8)  allows  the  exten¬ 
sion  of  Stokes’s  law  to  aerosol  size  below  0.01  Jim.  It  is  important  to  include  the  slip  correction  fac¬ 
tor  for  particles  whose  diameter  is  smaller  than  about  10  times  the  mean  free  path  of  gas  molecules. 


3.4  Gravitational  Settling 

Gravity  imposes  an  overall  downward  drift  on  particles,  which  contributes  an  enhanced  deposition 
flux  on  surfaces  with  an  upward  component  to  their  orientation.  The  gravity  force  exerted  on  a  parti¬ 
cle  is 


Fg  *”  ^  dp  (Pp  Pa  )§  ,  (9) 

where  g  is  the  acceleration  due  to  gravity  (9.81  m  s“2  at  the  Earth’s  surface),  and  pp  and  pa  are  particle 
and  air  densities,  respectively.  The  buoyancy  effect  usually  can  be  neglected  because  pa  is  commonly 
much  smaller  than  pp. 

The  terminal  settling  velocity  of  a  particle,  v/s,  is  established  when  the  gravity  force  is  balanced  by  the 
fluid  drag  force,  and  it  is  expressed  as27,28 


vts  = 


CcPpdp  8 
1 8|x 


(10) 


As  suggested  in  Eq.  (10),  the  particle  terminal  settling  velocity  increases  rapidly  with  particle  size 
since  it  is  proportional  to  the  square  of  particle  diameter  for  supermicron  particles  (where  Cc  =  1). 
Assuming  unit  particle  density,  for  instance,  the  settling  velocities  for  1  pm  and  10  pm  particles  are 
3.50  x  10'3  cm/s  and  3.05  x  10_l  cm/s,  respectively. 


3.5  Thermophoresis 

In  the  presence  of  a  temperature  gradient,  aerosol  particles  are  driven  from  the  high  to  low  tempera¬ 
ture  regions;  this  transport  process  is  known  as  thermophoresis.  The  phenomenon  was  first  reported 
in  the  19th  century,3^  and  its  quantitative  descriptions  were  published  by  Watson  in  193633  and  Zemik 
in  1957.34 

Thermophoretic  force  arises  from  asymmetrical  interactions  of  an  aerosol  with  the  surrounding  gas 
molecules  in  a  temperature  gradient,  where  gas  molecules  on  the  warm  side  bombard  the  particle  with 
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higher  average  momentum  than  those  on  the  cooler  side.  Acting  in  the  direction  of  decreasing  tem¬ 
perature,  thermophoretic  force  can  be  expressed  a$^ 


F,h=- 


37t^i2dpH  dT 
paT  d^’ 


(11) 


where  dT/dy  is  the  temperature  gradient,  and  H  is  the  thermophoretic  force  coefficient,  given  by 


H  = 


2.34 


ka/k, 


+ 


4.36  A/dj 


l  +  6.84A/dp  )( l  +  2ka/kp  +8.72A/dp 


(12) 


where  ka  and  kp  are  the  thermal  conductivities  of  air  and  the  particle  material,  respectively. 

When  the  thermophoretic  force  on  a  particle  is  balanced  with  fluid  drag,  the  thermophoretic  velocity 
is  obtained  as 


vth  =- 


CcvH  dT 
T  dy 


(13) 


Contrary  to  gravitational  settling  and  Brownian  diffusion,  which  are  strong  functions  of  particle  size, 
the  thermophoretic  velocity  is  nearly  independent  of  particle  size  for  particles  smaller  than  1  pm. 
This  occurs  because  both  the  thermophoretic  force  and  fluid  viscous  drag  have  approximately  the 
same  dependence  on  particle  diameter.  For  larger  particles,  the  thermophoretic  velocity  decreases 
with  increasing  particle  size  in  a  manner  that  depends  on  the  relative  thermal  conductivities  of  the 
particle  and  of  air. 


3.6  Electrostatic  Force 

A  charged  particle  migrates  in  an  electrical  field  due  to  the  Coulomb  force  Fc,  which  is  given  by 


Fc  =qE  =  neeE  , 


(14) 


where  q  is  the  charge  on  the  particle,  E  is  the  electric  field  strength,  ne  is  the  number  of  electrons  of 
deviation  from  the  electrically  neutral  state  (including  sign),  and  e  is  the  charge  of  a  single  electron  (- 
1.6  x  10~19  C). 

In  the  absence  of  an  electrical  field,  charged  aerosols  will  migrate  towards  or  away  from  a  conducting 
surface,  owing  to  the  image  force,  dielectric  force,  and  dipole-dipole  force,  although  these  forces  are 
much  weaker  than  the  Coulomb  force.  The  overall  electrostatic  force  on  a  charged  particle  can  be 
predicted  by*6 
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fe  =  qE  - 


(15) 


q2  |  qEdp3  37t£0dp6E2 
167i£0y2  1 6y  ^  1 28y 4 


where  £<>  is  the  permittivity  of  air  (8.86  x  10~12  C2  N  1  rrf 2),  and  y  is  the  normal  distance  from  a  sur¬ 
face.  The  terms  on  the  right-hand  side  of  Eq.  (15)  account  for  the  Coulomb  force,  image  force, 
dielectric  force,  and  dipole-dipole  force,  respectively.  When  an  electrical  field  is  present,  the  Cou¬ 
lomb  force  dominates,  and  both  the  dielectric  and  dipole-dipole  forces  are  negligible.  In  the  absence 
of  an  electric  field,  however,  the  only  electrostatic  force  responsible  for  particle  motion  is  the  image 
force,  as  suggested  in  Eq.  (15).  The  image  force,  which  only  occurs  near  a  conducting  surface,  is 
always  directed  toward  a  surface  and  may  dominate  over  Brownian  diffusion  and  turbulent  dispersion 
when  extremely  close  to  a  surface.  On  electrically  insulating  materials,  charges  may  accumulate  on 
the  surfaces,  giving  rise  to  electric  fields,  which,  in  turn,  affect  deposition  of  charged  particles. 

It  should  be  noted  that  aerosol  particles  acquire  or  lose  their  charges  through  random  collisions  with 
airborne  ions,  which  are  formed  by  ubiquitous  ionizing  radiation.  In  the  absence  of  an  electrical  field, 
the  processes  of  aerosol  acquiring  and  losing  charges  will  eventually  lead  to  an  equilibrium  charge 
state  called  Boltzmann  equilibrium.  The  maximum  number  of  charges  carried  by  a  particle  depends 
on  the  size  of  the  particles.  Table  2  shows  the  distribution  of  charges  carried  by  aerosols  at  Boltz¬ 
mann  equilibrium.  Take  10-pm  particles,  for  example.  At  equilibrium,  only  4.3%  of  the  particles  are 
electrically  neutral,  and  nearly  70%  of  the  aerosol  particles  carry  more  than  3  charges  on  a  single  par¬ 
ticle  (either  positive  or  negative). 


For  highly  charged  aerosols  in  an  electrical  field,  the  drift  velocity  resulting  from  electrostatic  forces 
on  particles  can  be  thousands  of  times  greater  than  that  from  the  gravitational  force. 


Table  2.  Distribution  of  Charge  on  Aerosols  at  Boltzmann  Equilibrium  (Reproduced  with 
permission  from  Hinds,  1999,  John  Wiley  &  Sons,  Inc.28) 


Particle 


Percentage  of  Particles  Carrying  the  Indicated  Number  of  Charges 


(jun) 

<  -3 

-3 

-2 

-1 

0 

+1 

+2 

+3 

>  +3 

0.01 

0.3 

99.3 

0.3 

0.02 

5.2 

89  6 

5.2 

0.05 

0.6 

19  3 

60.2 

19  3 

06 

0.1 

0.3 

4.4 

24.1 

42.6 

24.1 

4.4 

0.3 

0.2 

0.3 

2.3 

9.6 

22.6 

30.1 

22.6 

9.6 

2.3 

0.3 

0.5 

4.6 

6.8 

12.1 

17.0 

19.0 

17.0 

12.1 

6  8 

4.6 

1.0 

11.8 

8.1 

10.7 

12.7 

13.5 

12.7 

10.7 

8.1 

118 

2.0 

20.1 

7.4 

8.5 

9.3 

9.5 

9.3 

8.5 

7.4 

20.1 

50 

29.8 

5.4 

5.8 

6.0 

6.0 

6.0 

5.8 

54 

29  8 

10.0 

35.4 

4.0 

4.2 

4.2 

4.3 

4.2 

4.2 

40 

35.4 

12 


3.7  Turbophoresis 

Turbophoresis  is  a  phenomenon  in  which  the  net  migration  of  particles  occurs  from  regions  of  high  to 
low  turbulence  intensity,  i.e.,  toward  surfaces  where  the  turbulent  velocity  fluctuations  decrease  to 
zero.  Physically,  particles  in  regions  of  high  turbulence  intensity  acquire  sufficient  fluctuating  veloc¬ 
ity  components,  which  enable  them  to  drift  to  regions  where  the  turbulence  intensity  is  too  low  to 
send  them  back  with  sufficient  momentum  for  the  return  journey. 

Based  on  the  analogy  between  Brownian  motion  and  turbulent  diffusion,  as  well  as  energy  transfer 
from  the  fluid  to  the  particle,  the  mathematical  expression  for  turbophoretic  velocity  Vt  was  first  pro¬ 
posed  by  Caporaloni  et  ai  7  as 


V, 


^cPp^p  d(vy  ) 

18|i  dy 


(16) 


where  v'y  is  the  fluctuating  particle  velocity  normal  to  the  surface.  Reeks38  deduced  the  term  for  the 

turbophoretic  effect  and  arrived  at  the  same  expression  in  a  much  more  rigorous  approach  (a  special 
closure  of  Liouville  particle  equation  of  motion). 

Distinctly  different  from  turbulent  diffusion,  turbophoresis  is  attributed  to  the  interaction  between 
particle  inertia  and  the  inhomogeneous  turbulent  flow  field,  and  the  mass  transfer  takes  place  against 
the  gradient  in  turbulence  intensity,  even  in  the  absence  of  a  concentration  gradient.  The  effect  of 
turbophoresis  is  only  significant  for  particles  with  sufficiently  high  inertia,  and  is  typically  only  rele¬ 
vant  in  the  vicinity  of  a  surface  where  the  gradient  in  turbulence  intensity  is  high. 

For  near-surface  inertial  particles  under  highly  turbulent  flow  scenarios,  particle  transport  models 
accounting  for  turbophoresis  have  shown  good  agreements  with  experimental  measurements.39^*1 
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4.  Parameters  for  Particle  Deposition  Characterization 


To  characterize  the  rate  of  particle  deposition  in  an  enclosure,  deposition  velocity  (vd)  and  particle 
deposition  rate  (J5)  are  the  two  parameters  commonly  reported  in  the  literature.  Here,  their  defini¬ 
tions,  derivations,  and  physical  implications,  as  well  as  the  distinctions  between  (3  and  vdare 
discussed. 


4.1  Deposition  Velocity,  vd 

Particle  deposition,  a  mass  transfer  process  of  airborne  particles  onto  surfaces,  occurs  by  a  first-order, 
irreversible  mechanism.  Being  “first-order”  suggests  that  particle  deposition  flux  is  linearly  propor¬ 
tional  to  the  airborne  concentration,  and  being  “irreversible”  means  that  particles  simply  adhere  once 
they  collide  with  surfaces. 

The  aerosol  deposition  rate  can  be  parameterized  in  terms  of  a  mass  transfer  coefficient,  known  as  the 
deposition  velocity,  vd.  It  is  considered  as  the  proportionality  constant  between  the  net  aerosol  flux  J 
and  the  free-stream  airborne  concentration,  CM. 

J  r  i mass/ area  •  time  r  i length 

vd  =  — W - — - H— — •  (17) 

Coo  mass /volume  time 


The  aerosol  concentration,  ,  is  determined  at  a  position  sufficiently  far  away  from  the  surface,  i.e., 
core  zone  in  an  enclosure,  so  that  the  concentration  should  not  vary  greatly  with  positions. 

With  dimensions  of  length  per  time,  deposition  velocity  appears  to  represent  an  effective  velocity, 
which  incorporates  all  of  the  complexities  of  the  particle  deposition  process.  In  other  words,  deposi¬ 
tion  velocity,  as  an  aggregated  term,  comprises  every  aspect  of  particle  deposition  processes,  includ¬ 
ing  (1)  aerodynamic  transport  of  particles  by  turbulent  diffusion  from  the  well-mixed  core  region  to 
the  thin  boundary  layer  of  air  adjacent  to  the  surface;  and  (2)  mass  transfer  of  particles  across  the 
boundary  layer,  and  subsequent  uptake  by  the  surface  through  a  variety  of  particle  deposition 
mechanisms. 

The  magnitude  of  the  deposition  velocity  depends  on  factors  that  govern  particle  transport:  particle 
size,  the  near-surface  air  turbulence,  surface  characteristics  including  orientation  and  texture  (smooth 
vs.  rough),  air-surface  temperature  difference,  and  the  presence  of  an  electrical  field  near  the  surface. 
A  literature  review  of  the  existing  experimental  findings  concerning  the  influence  of  these  factors  on 
particle  deposition  is  presented  in  Section  6. 
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4.2  Particle  Deposition  Rate,  p 

Another  common  parameter  used  to  quantitatively  characterize  the  particle  deposition  rate  (3  in  an 
enclosure  is  defined  as:20 


dC(dp,l) 

dt 


=  -pc. 


(18) 


where  C(dp,t)  is  the  aerosol  number  or  mass  concentration  as  a  function  of  time,  t,  in  the  core  region 
of  the  enclosure,  and  (3  is  the  particle  deposition  rate  with  the  unit  of  time"1  (e.g.,  s"1).  Particle 
deposition  rate  is  also  commonly  referred  to  as  the  aerosol  removal  rate  or  wall  loss  rate  because  the 
deposition  process  removes  particles  from  the  gas  phase  onto  the  available  wall  surfaces. 

Equation  (18)  is  written  in  the  form  of  mass  balance  of  particles  suspended  in  air  within  an  enclosure, 
with  the  following  conditions:  (1)  no  or  negligible  sources  for  particle  generation  (such  as  particle 
influx  in  the  incoming  air  flow  and  nucleation);  (2)  no  particle  size  change  due  to  coagulation,*  con¬ 
densation,  and  evaporation;  and  (3)  deposition  onto  surfaces  is  the  only  particle  loss  mechanism. 

Note  that  the  particle  deposition  rate  is  expected  to  be  a  function  of  particle  size,  as  is  the  case  for 
particle  deposition  velocity.  Also  note  that  Eq.  ( 1 8)  only  accounts  for  particle  loss  attributable  to 
deposition.  In  an  enclosure,  there  are  typically  other  processes  occurring  in  parallel,  such  as  air 
exchange,  that  influence  the  particle  concentration  and  should  be  taken  into  account.6 

Experimentally,  the  aerosol  concentration  decays  with  time  as  expressed  by  Eq.  (18)  after  initial 
mixing  of  aerosols  in  an  air-tight  enclosure.  Mathematically,  the  rearrangement  of  Eq.  ( 1 8)  gives 

c,  -C0e-*\  (19) 

where  C0  is  the  initial  aerosol  concentration  in  the  enclosure,  Q  is  the  aerosol  concentration  at  time  t. 
When  the  air  is  well-mixed,  plotting  InC  as  a  function  of  time  gives  rise  to  a  linear  plot,  with  the 
negative  of  the  slope  yielding  the  particle  deposition  rate,  (3.  Figure  2  represents  an  example  of  parti¬ 
cle  concentration  decay  as  a  function  of  time  and  the  associated  deposition  rate.  Similar  to  particle 
deposition  velocity  vd,  (3  incorporates  all  the  deposition  processes  that  remove  aerosols  from  being 
suspended  in  the  enclosure. 


4.3  Comparison  of  vd  and  (3 

The  major  distinction  between  vd  and  (3  is  that  vd  depends  on  surface  orientations  with  respect  to 
gravity  (a  “local”  parameter),  while  (3  represents  an  average  term  over  all  available  surfaces  for  depo¬ 
sition  within  an  entire  enclosure  (a  “global”  parameter).  For  instance,  in  the  context  of  a 


Particle  coagulation,  which  contributes  to  aerosol  size  distribution  change,  occurs  when  the  aerosol  concentration  is 
sufficiently  high  When  that  occurs,  coagulation  may  be  an  important  mechanism  to  incorporate  in  Eq.  (18),  as  seen  in 
Okuyama  et  alP 
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Figure  2.  An  illustrative  example  of  aerosol  concentration  decay 
measurements  of  various  particle  diameters  as  a  function 
of  time  in  a  165-L  cylindrical  chamber  (reprinted  with 
permission  from  Chen  et  alA1).  The  particle  deposition 
rates,  obtained  from  the  values  of  the  slope,  were  0.058, 

0.14,  0.088,  0.28,  0.44,  0.71  for  particle  diameters  of 
0.305,  0.073,  0.91,  1.74,  2.02,  and  2.99  pm,  respectively. 

rectangular  enclosure,  the  magnitude  of  vd  for  supermicron  particles  varies  strongly  with  surface  ori¬ 
entation,  i.e.,  vdtfloor  >  vd  wan  >  vdtCci]ing,  owing  to  the  dominant  contribution  of  gravitational  settling  over 
diffusion  and  other  mechanisms. 

The  value  of  (3  can  be  related  to  the  values  of  vd  for  all  enclosure  surfaces.  By  material  balance,  the 
deposition  rate,  P,  is  the  surface-area-averaged  vd  multiplied  by  the  surface-to-volume  ratio  (S/V)  of 
the  enclosure: 


^  ~  Jvd  (s)ds  =  ■ -C  v  d  =  -PC  ,  (20) 

where  vd(s)  is  the  particle  deposition  velocity  onto  the  surface  s.  The  integration  is  to  be  carried  out 
over  all  the  enclosure  interior  surfaces  with  the  total  surface  area  S,  and  vd  is  the  area-weighted  mean 
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deposition  velocity.  It  can  be  clearly  seen  that  (3,  unlike  vd,  eliminates  the  explicit  spatial  dependence 
of  particle  deposition. 
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5.  Methods  for  Measurement  of  Particle  Deposition 


Owing  to  its  importance  for  numerous  engineering  systems,  the  processes  of  particle  deposition  on 
surfaces  have  been  widely  studied  experimentally.  Here,  the  general  approach  for  measuring  (3  and  vd 
with  respect  to  an  enclosed  system  is  presented  with  the  aim  of  helping  to  gain  insights  into  the 
experimental  findings  that  will  be  presented  in  Sections  6  and  7,  respectively. 


5.1  Measuring  Particle  Deposition  Rate,  p 

Experimentally,  the  particle  deposition  rate  P  as  a  function  of  particle  diameter  is  inferred  from  aero¬ 
sol  concentration  decay  over  time,  following  Eq.  (18)  and  as  shown  in  Figure  2.  To  do  so,  monodis- 
perse  aerosols  are  usually  generated  in  elevated  concentrations  so  that  sufficient  time  is  available  to 
take  multiple  data  points  on  aerosol  concentrations  during  the  course  of  deposition  experiments. 
When  the  test  aerosols  are  polydisperse,  the  aerosol  measuring  instruments  will  need  to  have  the  par¬ 
ticle  size-resolved  capability  in  order  to  determine  the  concentration  decay  rate  as  a  function  of  parti¬ 
cle  diameter.  The  operation  principles  for  some  common  particle  sizing  instruments,  e.g.,  optical 
particle  counter,  time-of-flight  aerosol  spectrometer,  electrical  aerosol  mobility  analyzer,  as  well  as 
aerosol  generation  methods,  can  be  found  in  references  [28,43,44]. 

Note  that  the  degree  of  air  turbulence  plays  an  important  role  in  influencing  the  rate  of  particle  depo¬ 
sition.  Therefore,  the  stirring  of  air  within  the  experimental  chamber,  by  either  impeller  mixing  or 
flushing,  should  be  controlled  consistently  throughout  the  experiment. 

The  scheme  of  determining  particle  deposition  rate  by  monitoring  airborne  particle  concentration 
decay  over  time  does  not  involve  direct  measurements  of  deposited  particles  on  surfaces,  thus  the 
measured  deposition  rate  is  “inferred.”  This  type  of  experiment  is  relatively  simple  to  carry  out,  and 
data  analysis  is  straightforward.  However,  the  drawback  of  this  indirect  approach  is  the  inability  to 
differentiate  the  contribution  of  particle  deposition  on  surfaces  with  respect  to  orientations  and 
locations. 


5.2  Measuring  Particle  Deposition  Velocity,  vd 

As  suggested  in  Eq.  (17),  particle  deposition  velocity  is  determined  by  normalizing  the  deposition 
flux  with  the  average  aerosol  concentration  over  the  duration  of  particle  deposition.  Distinct  from 
measuring  the  deposition  rate  constant,  P,  evaluating  vd  requires  a  direct  measurement  of  particle 
deposition,  which  involves  the  simultaneous  determination  of  particle  deposition  flux  and  airborne 
particle  concentration.  In  other  words,  the  deposited  particle  mass  will  need  to  be  recovered  from  the 
surface  (or  otherwise  detected  and  quantified),  and  the  aerosol  concentration  outside  the  boundary 
layer  during  the  course  of  deposition  will  also  need  to  be  determined.  For  the  typical  experiment  of 
this  type,  monodisperse  aerosols  need  to  be  generated  and  injected  into  the  experimental  chamber  to 
provide  the  size-specific  deposition  flux.  To  obtain  the  deposition  velocity  as  a  function  of  particle 
diameter,  the  same  experimental  procedures  are  repeated  for  other  particle  sizes  of  interest.  The 
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advantage  of  this  direct  measurement,  albeit  time-consuming  and  labor-intensive,  is  that  spatially 
resolved  particle  deposition  can  be  determined.  A  key  challenge  of  this  approach,  however,  is  the 
need  to  accurately  determine  trace  amounts  of  particles  deposited  on  surfaces.  The  earliest  method 
was  performed  by  microscopically  counting  the  particles  on  a  small  area  of  a  test  surface.45  Despite 
the  recent  rapid  advancement  of  optical  and  digital  imaging  techniques,  it  remains  a  tedious  task  to 
perform  the  counting,  considering  that  only  a  very  small  surface  area  is  available  per  microscopic 
frame,  which,  in  turn,  limits  the  possibility  of  sampling  larger  surface  areas. 

To  quantitatively  determine  deposited  particle  mass  on  surfaces,  aerosols  can  be  labeled  to  facilitate 
subsequent  surface  analysis.  Two  types  of  tracer  aerosol  detection  techniques  have  been  employed 
for  such  purposes:  (1)  fluorescence  spectroscopy,  and  (2)  neutron  activation  analysis  (NAA)  or  pro¬ 
ton-induced  X-ray  emission  (PIXE).  With  respect  to  the  former  method,  monodisperse  aerosols  are 
generated  with  fluorescent  materials.  After  deposition,  the  fluorescent  particles  on  the  test  surface  are 
extracted  with  a  solution  of  known  volume.  The  resulting  fluorescence  intensity  of  the  solution  is 
proportional  to  the  collected  particle  mass  from  the  surface.46-48  In  order  to  accumulate  sufficient 
fluorescent  particle  mass  for  analysis,  the  experimental  deposition  time,  particularly  for  submicron 
particles,  could  easily  exceed  100  h.46 

Contrary  to  the  chemistry-based  fluorescence  technique,  NAA  or  PIXE  relies  on  the  physics  of  atoms 
in  the  aerosol  materials.  In  NAA,  the  energy  of  the  induced  radioactivity  through  neutron  bombard¬ 
ment  is  characteristic  of  a  particular  element,  while  PIXE  involves  the  measurement  of  the  character¬ 
istic  X-ray  emission  via  high-energy  proton  excitations  of  the  elements  in  the  sample.  Methods  using 
NAA  and  PIXE  are  semi-invasive;  i.e.,  only  the  particle-bearing  surface  needs  to  be  removed  for  the 
analysis  from  the  test  chamber,  and  removal  of  deposited  particles  from  the  particle-bearing  surface  is 
not  required.  To  minimize  interfering  effects  from  various  elements,  rare  earth  elements  such  as  dys¬ 
prosium  (164Dy)  and  indium  (ll5In)  can  be  incorporated  into  the  aerosols  for  the  purpose  of  subsequent 
NAA  analysis  49-51  Owing  to  the  extremely  high  sensitivity,  for  instance,  a  l64Dy  mass  of  the  order  of 
10  10  g  can  be  easily  detected;  hence,  the  experimental  deposition  time  can  be  significantly  reduced 
(e.g.,  15-20  min).50  The  drawback  of  this  technique,  however,  is  that  specialized  facilities,  such  as  a 
nuclear  reactor  for  neutron  activation,  is  required. 


20 


6.  Review  of  Experimental  Studies 


A  significant  body  of  experiments  has  been  devoted  to  exploring  particle  deposition  in  enclosures 
since  the  late  1940s  due  to  the  relevance  to  material  deterioration  and  the  implications  to  human 
health  risk,  etc.  Table  A1  in  the  Appendix  highlights  a  chronological  summary  of  the  laboratory 
investigations  pertaining  to  particle  deposition  in  an  enclosed  environment. 

The  processes  of  particle  deposition  onto  surfaces  are  complex,  multifaceted  phenomena  since  they 
vary  strongly  depending  on  the  characteristics  of  particles  and  airflow  patterns,  as  well  as  surface 
properties.  In  brief,  the  processes  of  particle  deposition  consist  of  two  mechanisms  in  series:  (1)  par¬ 
ticle  transport  from  the  core  region  in  the  enclosure  to  the  boundary  layer  adjacent  to  the  surface,  and 
(2)  subsequent  deposition  onto  the  surface.  Consequently,  both  properties  of  airflow  and  the  surface 
play  a  crucial  role  in  governing  the  rate  of  particle  deposition  in  steps  (1)  and  (2),  respectively.  The 
degree  of  airflow  turbulence  controls  how  rapidly  particles  migrate  to  the  proximity  of  a  surface, 
while  the  surface  characteristics,  such  as  orientations  and  roughness,  determine  how  readily  particles 
interact  with  the  surface.  As  will  be  seen  soon,  particle  size  is  the  most  important  parameter  govern¬ 
ing  the  motions  of  particles  in  a  gas  phase  associated  with  the  transport  properties. 

Here,  the  important  physical  factors  that  influence  particle  deposition  onto  enclosure  surfaces  are 
summarized  from  the  existing  literature. 


6.1  Particle  Characteristics  Affecting  Particle  Deposition 

Particle  size  is  key  in  the  determination  of  particle  deposition.  The  dimension  of  particles  in  the  gas 
phase  is  commonly  characterized  in  terms  of  aerodynamic  diameter,  which  refers  to  an  equivalent 
diameter  of  a  unit  density  spherical  particle  with  the  same  terminal  settling  velocity  as  the  particle 
being  measured/  Therefore,  the  concept  of  aerodynamic  diameter  has  taken  particle  shape  and  den¬ 
sity  into  account,  and  it  has  successfully  captured  the  aerodynamic  behavior  of  airborne  particles  in 
many  systems  of  interest/  The  notion  of  aerodynamic  diameter  has  been  commonly  used  in  the  aero¬ 
sol  literature,  including  this  review. 

Figure  3  illustrates  the  characteristic  V-shaped  curves  of  experimentally  determined  particle  deposi¬ 
tion  rates  as  a  function  of  particle  diameter.  The  large  variability  of  the  data,  which  can  be  up  to  2 
orders  of  magnitude,  reflects  the  fact  that  other  factors,  such  as  airflow  turbulence  and  surface  char¬ 
acteristics,  also  play  a  significant  role  in  influencing  the  extent  of  particle  deposition  in  enclosures. 

As  mentioned  earlier,  when  particles  are  not  electrically  charged,  the  rate  of  particle  deposition  is 
governed  by  gravitational  settling  and  Brownian  motion  of  particles,  as  well  as  the  turbulent  motion 
of  the  fluid  in  which  particles  are  entrained.  For  particles  smaller  than  0. 1  pm,  Brownian  and  turbu- 


The  behavior  of  nonspherical  particles  in  shear  flows  might  not  be  properly  characterized  using  the  concept  of 
aerodynamic  diameter  owing  to  their  complex  rotational  and  translational  motions  not  accounted  for  in  spherical  particles. 
Thus,  the  prediction  and  measurement  of  deposition  using  aerodynamic  diameters  for  nonspherical  particles  under  this 
scenario  would  be  an  average  outcome  contributing  from  their  stochastic  behaviors  in  the  turbulent  flow. 
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Particle  diameter,  pm 

Figure  3.  Comparison  of  experimentally  obtained  particle  deposition  rates  as  a  function  of  panicle 
diameter.  The  values  of  particle  deposition  rate  p  exhibit  distinct  minima  for  0.1-1  pm 
particles,  as  indicated  in  the  V-shaped  curves. 


lent  diffusion  control  particle  deposition.  On  the  other  hand,  gravitational  settling  becomes  dominant 
in  removing  particles  from  the  gas  phase  with  the  increase  of  particle  size,  e.g.,  >1  pm.  As  a  result,  a 
minimum  of  particle  deposition  rate  occurs  for  particles  of  0.1-1  pm,  which  are  generally  referred  to 
“accumulation  mode”  particles,  when  neither  of  the  mechanism  works  effectively  to  cause  particle 
deposition  (Figure  3). 


With  respect  to  particle  composition,  experiments  have  shown  that  the  material  of  particles  makes  no 
difference  in  the  measured  particle  deposition  rates. Also,  as  seen  in  Table  Al,  various  aerosol 
materials  have  been  employed  for  particle  deposition  experiments,  and  the  choice  of  aerosol  compo¬ 
sitions  does  not  alter  the  outcome. 


6.2  Airflow  Characteristics  Affecting  Particle  Deposition 

Carried  by  large-scale  eddy  currents,  airborne  particles  are  brought  from  one  point  to  another  by  tur¬ 
bulent  diffusion.  Eddy  currents  can  be  induced  by  natural  convection  (e.g.,  temperature  gradient), 
and  forced  convection  (e.g.,  mechanical  stirring  and  fluid  flushing). 


6.2.1  Forced  Convection 

Early  experimenters  observed  that  increasingly  turbulent  air  motion  in  an  enclosed  space  was  related 
to  enhanced  particle  deposition  onto  surfaces,  and  the  empirical  parameters  describing  such  phenom¬ 
ena  were  obtained  as  a  function  of  the  air  stirring  intensity.52  Corner  and  Pendlebury53  later  explained 
these  empirical  observations  based  on  the  theoretical  grounds,  and  concluded  that  air  turbulence  was 
one  important  factor  that  influenced  the  rate  of  particle  deposition. 
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Subsequent  experiments  help  to  shed  light  on  the  dependence  of  particle  deposition  rates  on  the  air 
turbulence  level  in  the  enclosures.  Figure  4  represents  four  laboratory  results  of  particle  deposition 
rates  with  respect  to  different  air  mixing  scenarios.  It  can  be  seen  that  enhanced  particle  deposition  is 
attributed  to  increased  air  turbulence,  which  was  produced  by  either  fluid  flushing  (Nomura  el  a/.54) 
or  fan  stirring  (Okuyama  el  al .23,  Shimada  el  a/.55,  and  Cheng56).  The  enhancement  of  particle  deposi¬ 
tion  onto  enclosure  surfaces  results  from  more  effective  aerodynamic  mass  transfer  at  higher  airflow 
turbulence,  which  brings  aerosol  particles  more  rapidly  toward  a  surface. 


In  addition,  as  the  turbulent  fluctuation  of  air  motion  becomes  significant,  particle  inertia  may  con 
tribute  to  particle  deposition,  and  inertial  transport  of  particles  through  the  boundary  layer  could  be 
potentially  enhanced  by  near-surface  turbulent  bursts.41  Inertially  induced  turbulent  deposition  could 
be  particularly  important  for  large  particles. 
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Figure  4.  Dependence  of  experimentally  measured  particle  deposition  rates  on  the  degree  of  airflow 
turbulence,  as  characterized  by  the  fan  stirring  rate  (in  rpm)  and  the  air-exchange  rate 
(incoming  airflow  rate  divided  by  the  enclosure  volume,  with  the  unit  of  h_1).  The  interior 
surface  area-to- volume  ratio,  S/V,  is  indicated  for  each  test  chamber.  Within  the  same 
experimental  system,  enhanced  particle  deposition  is  consistently  observed  at  higher  air 
turbulence  levels. 


23 


6.2.2  Natural  Convection 

Natural  convection  arises  from  temperature  differences  among  air  parcels,  or  heat  transfer  at  surfaces 
(i.e.,  surface-to-air  temperature  difference).  In  the  absence  of  forced  convection,  natural  convection 
becomes  the  only  means  of  air  mixing  inside  enclosed  spaces.  Cheng56  conducted  particle  deposition 
experiments  under  isothermal  and  still  air  conditions  with  no  apparent  air  movements  detected  in  the 
chamber,  but  still  found  the  deposition  data  well  described  by  the  homogeneous  turbulence  model, 
assuming  the  core  region  was  well-mixed  (see  Subsection  7.1  for  more  discussion). 

Figure  5  shows  a  comparison  of  particle  deposition  rates  as  a  function  of  particle  size  under  natural 
convection  conditions  from  different  experiments.27  42  5fv  ss  First,  notice  that  particle  size  remains  the 
predominant  parameter  governing  the  deposition  rates,  as  indicated  by  the  distinct  V-shaped  curves. 
Secondly,  the  scatter  of  the  data  in  these  experiments  may  be  attributed  to  (1)  the  S/V  ratio,  as  higher 
surface  area  normalized  by  volume  translates  to  more  surfaces  available  for  particles  to  deposit  on, 
and  (2)  temperature  gradient  within  the  enclosure,  as  it  dictates  the  extent  of  convective  mixing. 

Chen  et  alA  and  Chengs*  specifically  documented  the  temperature  measurements  inside  their  test 
chambers,  ensuring  that  the  experimental  systems  were  either  isothermal  or  otherwise  reported.  The 
temperature  data  from  the  other  deposition  measurements  are  either  insufficient  or  unavailable  to 
draw  any  useful  information. 


6.3  Surface  Characteristics  Affecting  Particle  Deposition 

After  airborne  particles  enter  a  boundary  layer,  they  may  be  transported  to  the  surface  by  means  of  a 
variety  of  deposition  mechanisms,  as  described  in  Section  3.  It  remains  a  good  approximation  and 
has  been  demonstrated  experimentally  that  adhesion  of  micrometer-sized  (or  smaller)  particles  is 


Particle  diameter,  pm 

Figure  5.  Comparison  of  experimentally  measured  particle  deposition  rates  as  a  function  of  particle 
size  under  natural  convection  conditions.  The  interior  surface  area  to  volume  ratio,  S/V,  is 
indicated  for  each  experimental  chamber.  AT  in  Chen  et  al.A 2  refers  to  the  temperature  dif¬ 
ference  between  the  top  and  bottom  walls  (e.g.,  Tb-  Tt  =  10°C)  of  the  test  chamber. 


*  The  detection  limit  in  these  experiments  was  about  1  cm/s. 
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complete  and  irreversible  once  they  come  into  contact  with  a  surface.  The  Van  der  Waals  force  is 
the  predominant  adhesion  force  between  an  aerosol  particle  and  any  surface  (including  another  parti 
cle),  although  electrostatic  and  capillary  forces  may  be  important  as  wellr 

Surface  characteristics  play  a  role  influencing  particle  deposition  as  they  often  affect  particle-surface 
interactions  within  the  boundary  layer.  The  dependence  of  particle  deposition  on  surface  characteris¬ 
tics  has  been  explored  experimentally  by  manipulating  various  factors  such  as  surface  orientations 
and  roughness.  Below,  an  overview  is  provided  concerning  experimental  deposition  studies  with 
respect  to  various  surface  characteristics. 


6.3.1  Surface  Orientation:  Horizontal  vs.  Vertical 

The  effect  of  gravity  contributes  to  the  major  difference  in  terms  of  particle  deposition  flux  onto  sur¬ 
faces  of  various  orientations,  e.g.,  a  horizontal  upward-facing  surface  (floor)  as  opposed  to  a  vertical 
(wall)  or  a  horizontal  downward-facing  surface  (ceiling).  Briefly,  the  mechanism  of  gravitational 
settling  is  primarily  responsible  for  particle  deposition  onto  the  floor,  while  turbulent  and  Brownian 
diffusions  are  dominant  deposition  processes  to  walls  and  ceiling.  Particle  flux  with  respect  to  sur¬ 
faces  of  various  orientations  can  only  be  studied  by  means  of  spatially  resolved  deposition  velocity 
measurements  as  described  in  Subsection  5.2. 

Byrne  et  al.49  employed  the  NAA  technique  to  recover  the  particle  mass  deposited  on  various  surfaces 
under  forced  convection  in  their  cubic  test  chamber  (2x2x2  nr ).  They  reported  that  proportion¬ 
ately  more  of  the  overall  particle  deposition  flux  was  found  on  the  floor  (and  thus  less  on  walls)  as 
particle  size  increases,  as  summarized  in  Table  3. 

Thatcher  et  al ,46  used  fluorescent  tracer  particles  to  measure  particle  deposition  velocity  with  respect 
to  different  surfaces  under  natural  convection  flow  conditions  in  their  cubic  enclosure  (1.2  x  1 .2  x  1.2 
m^).  In  addition  to  examining  the  effect  of  different  surface  orientations,  the  interior  surface  tem¬ 
peratures  of  the  chamber  were  independently  controlled  at  fixed  surface-to-air  temperatures  to  inves¬ 
tigate  the  thermophoretic  effects  on  particle  deposition  (see  Subsection  6.3.3  for  more  discussion).  In 
brief,  the  experimental  chamber  was  heated  on  the  floor  and  one  vertical  wall,  and  cooled  on  the 
ceiling  and  the  opposite  wall.  As  long  as  the  surface-to-air  temperature  difference  is  kept  identical, 
the  thermophoretic  effect  is  considered  to  be  very  similar,  allowing  comparisons  to  be  made  con¬ 
cerning  particle  deposition  on  surfaces  of  different  orientations. 


Table  3.  Relative  Contributions  of  Particle  Deposition  Flux  Experimentally  Determined  for  the 
Horizontal  and  Vertical  Surfaces  as  a  Function  of  Particle  Diameter  in  Byrne  et  al4" 


Particle  size 
Oun) 

Flux  on  horizontal  surface 
(floor) 

Flux  on  vertical  surface 
(one  wall) 

Average  deposition  velocity* 
(m/s) 

0.7 

57% 

9% 

4.1x10  s 

2.5 

68% 

8% 

6.2x10  s 

4.5 

72% 

7% 

1.1x10"* 

5.4 

80% 

5% 

2.0x1 0* 

*  accounted  for  particle  deposition  to  all  chamber  surfaces 
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Figure  6  presents  a  comparison  of  particle  deposition  velocity  experimentally  obtained  from  vertical 
and  horizontal  surfaces  at  the  same  temperature  [+1 .5K  in  (a)  and  -1.5K  in  (b)  relative  to  the  air  tem¬ 
perature].  One  of  the  distinctive  features  in  Figure  6(a)  is  that  particle  deposition  velocity  onto  the 
horizontal  upward-facing  surface  (i.e.,  floor)  increases  as  particle  size  increases,  in  agreement  with 
theoretical  predictions  and  experimental  results  in  Byrne  et  alf  Also  note  in  Figure  6  (a)  that  the 
measured  deposition  velocities  onto  vertical  surfaces  (i.e.,  wall),  were  increasingly  lower  than  those 
on  the  floor  by  1  to  3  orders  of  magnitude  with  the  increase  of  particle  size.  Due  to  the  diminishing 
contribution  from  gravitational  settling,  the  deposition  velocity  to  the  horizontal  downward-facing 
surface  (i.e.,  ceiling)  also  decreases  for  increasing  particle  size,  as  shown  in  Figure  6(b). 


Particle  diameter,  jim 


Figure  6. 


Particle  diameter,  nm 

Comparison  of  experimentally  measured  deposition  velocity  for  surfaces  of 
various  orientations  in  a  cubic  chamber  under  natural  convection.  The  surface 
temperatures  in  (a)  and  (b)  were  kept  at  1 .5K  higher  and  1 ,5K  lower  than  that 
of  air,  respectively.  The  values  of  deposition  velocity  plotted  here  were 
obtained  from  Table  2  in  Thatcher  et  al4f  The  symbol  of  *  in  (b)  denotes  that 
the  deposited  particle  mass  on  the  ceiling  was  under  the  detection  limit. 
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6.3.2  Surface  Texture:  Smooth  vs.  Rough 

Surface  texture  complicates  the  process  of  particle  deposition.  The  microscale  roughness  elements 
can  alter  the  near-surface  turbulent  structures  and  reduce  the  boundary  layer  thickness,  which,  in  turn, 
influences  particle  deposition.  Roughness  can  also  influence  diffusive  deposition  even  when  the 
fluid-mechanical  properties  are  not  disturbed.  The  roughness  elements  can  extend  into  and  across  the 
particle  concentration  boundary  layer  (which  can  be  much  thinner  than  the  viscous  sublayer),  and  this 
can  expose  the  elements  to  higher  particle  concentration,  thus  enhancing  the  deposition  rate.  In  addi¬ 
tion,  particle  deposition  may  be  further  enhanced  by  more  sites  within  the  roughness  elements  avail¬ 
able  for  impaction  provided  that  particle  inertia  is  sufficient.  Experimental  data  regarding  the  effect 
of  surface  roughness  on  particle  deposition  within  an  enclosure  are  sparse.  Nonetheless,  a  large  col¬ 
lection  of  literature  exists  for  particle  deposition  onto  smooth  and  rough  pipe  walls,  allowing  some 
clues  to  be  inferred.  Sippola 15  provides  an  excellent  review  of  this  issue. 


Based  on  some  assumptions  and  the  measurements  from  a  stirred  2.6-L  cylindrical  chamber,  Shimada 
et  tf/.6u  proposed  a  semi-empirical  model  to  explain  the  dependence  of  surface  roughness  height  on 
particle  Brownian  and  turbulent  diffusive  deposition,  and  to  estimate  the  rate  of  particle  deposition. 
Figure  7  illustrates  the  influence  of  surface  roughness  level  on  the  experimentally  measured  deposi¬ 
tion  rates  for  0.0 1-0.2  |im  particles  under  the  same  turbulent  mixing  conditions  (500  rpm).  They 
concluded  that  as  the  surface  roughness  becomes  significant,  particle  deposition  tends  to  be  influ¬ 
enced  by  the  turbulent  diffusion  very  close  to  the  surface,  leading  to  the  enhancement  of  particle 
deposition.  Later,  Shimada  et  al.62  further  refined  the  model,  and  reported  that  the  experimentally 
determined  deposition  rates  agreed  well  with  those  reproduced  in  model  calculations  when  the  two- 
dimensional  configuration  of  surface  roughness  was  taken  into  account  in  calculating  particle  con¬ 
centration  above  a  rough  surface.  Since  their  modeling  approach  focused  on  diffusive  deposition 
ofO.0 1-0.2  pm  particles,  the  mechanisms  of  gravitational  settling  and  inertial  impaction  were  not 
included. 


0.001  0.01  01  1  10 
Particle  diameter,  gm 

Figure  7.  The  dependence  of  experimentally  measured  particle  deposition  rates  on  the  surface  rough¬ 
ness.  The  length  scale  of  1 1 .2,  1 19.2,  and  204.8  pm  in  Shimada  et  al.60  refers  to  the  aver¬ 
age  height  of  roughness  of  the  sandpapers  placed  on  the  interior  surface  of  the  enclosed 
tank.  The  S/V  ratios  of  the  experimental  chamber  in  Shimada  et  and  Abadie  et  al.6]  are 
47  and  10  m_I,  respectively. 


*  Shimada  et  al.  [60],  1 1 .2  jam 
o  Shimada  et  al.  [60],  119.2  [am 
o  Shimada  et  al.  [60],  204.8  urn 

*  Abadie  et  al.  [61  ],  linoluem 

£  Abadie  et  al.  [61],  smooth  wall  paper 
a  Abadie  et  al.  [61],  rough  wall  paper 

*  Abadie  et  al.  [61  ],  carpet 
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Thatcher  and  Nazaroff47  studied  the  influence  of  different  surface  roughness  on  particle  deposition 
under  natural  convection  in  a  1 .8  cubic  enclosure.  They  found  that  deposition  of  small  particles 
(~0.2  pm)  was  relatively  insensitive  to  surface  textures,  but  more  pronounced  effects  were  observed 
with  increasing  particle  size.  In  their  experiments,  the  measured  deposition  onto  the  roughest  verti¬ 
cal*  surface  for  1 .3  pm  particles  was  5  times  greater  than  deposition  to  a  vertical  smooth  surface. 
They  also  noted  that  the  surface  roughness  had  more  substantial  effects  for  vertical  surfaces  than 
horizontal,  and  for  warm  surfaces  than  cool. 

Various  wall  treatments  on  all  interior  surfaces  were  furnished  in  a  cubic  chamber  (0.6  x  0.6  x  0.6 
rn)  to  measure  deposition  rates  in  Abadie  et  a! .61  for  particles  of  0.7,  1 .0,  and  5.0  pm  under  fan  mix¬ 
ing  conditions.  As  shown  in  Figure  7,  the  measured  particle  deposition  rates  were  found  to  increase 
with  the  following  order  of  increasing  surface  roughness  level:  linoleum,  smooth  wall  paper,  rough 
wall  paper,  and  carpet. 

Lai  et  al.6}  conducted  deposition  experiments  in  an  8-m*  test  chamber  under  different  fan  speeds 
using  monodisperse  particles  from  0.7  to  5.4  pm.  They  observed  that,  at  the  highest  fan  speed,  the 
average  ratio  of  particle  deposition  flux  onto  a  vertical  rough  surface  over  smooth  surface  increased 
from  1 .05  for  0.7  pm  particles  to  1.1,  1.6,  and  2.4  for  particle  sizes  2.5,  4.3,  and  5.4  pm,  respectively. 
Since  gravity  does  not  play  a  role  in  particle  deposition  onto  vertical  walls,  surface  roughness  and 
perhaps  particle  inertia  for  larger  particles,  are  believed  to  contribute  to  the  enhanced  particle  deposi¬ 
tion  in  the  experiments. 

Thatcher  et  a/.64  reported  on  the  influence  of  particle  deposition  rates  upon  different  furnishing  levels 
(bare,  carpeted,  and  fully  furnished  with  chairs,  curtains,  etc.)  in  a  room-sized  setting  (2.2  x  2.7  x  2.4 
m3)  for  particle  sizes  from  0.5  to  10  pm.  The  use  of  carpeting  and  furniture  can  be  considered  to 
increase  the  average  roughness  of  surfaces  (although  the  airflow  pattern  will  be  different  for  the  fur¬ 
nished  scenario  due  to  additional  obstructions  inside  the  room;  also  increased  surface  area  to  volume 
ratio).  Across  all  particle  sizes  measured,  the  measured  deposition  rates  were  observed  to  have  the 
consistently  descending  trend  with  respect  to  the  furnishing  or  roughness  level:  furnished  >  carpeted 
>  bare. 

Lai  and  Nazaroff48  placed  glass  plates  and  sandpapers  of  different  grades  on  the  vertical  walls  of  a 
1 .8-m^  cubic  chamber  to  measure  the  deposition  velocities  of  supermicron  particles  under  forced  con¬ 
vection  conditions.  They  reported  that  the  experimentally  determined  deposition  velocity  appeared  to 
reach  a  fairly  steady  value  for  particles  larger  than  7  pm.  In  addition,  the  deposition  velocity  was 
observed  to  increase  with  increasing  roughness  grade,  albeit  the  increments  with  increasing  surface 
roughness  were  not  as  significant  as  the  influence  of  particle  size  on  deposition. 


6.3.3  Surface  Temperature:  Warm  vs.  Cold 

As  discussed  in  Subsection  3.5,  thermophoretic  forces  on  airborne  particles  can  be  induced  by  tem¬ 
perature  differences  between  a  surface  and  air.  In  the  presence  of  a  temperature  gradient,  particles 
always  move  toward  the  direction  of  lower  temperature.  As  a  consequence,  particles  tend  to  prefer¬ 
entially  deposit  onto  a  cold  surface  over  a  warm  one  due  to  thermophoresis. 


Vertical  surfaces  are  studied  to  exclude  the  direct  influence  of  gravity  on  deposition 
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Numerous  experi  menial  studies  have  examined  the  effectiveness  of  exploiting  thermophoresis  to  mini¬ 
mize  particle  deposition  on  silicon  wafers  for  microcontamination  control.1  1  However,  experimental 
data  for  particle  deposition  on  cold  and  warm  surfaces  with  respect  to  an  enclosed  volume  are  sparse. 
The  experiments  performed  by  Thatcher  et  al .46  provided  valuable  insights  into  particle  deposition 
under  the  influence  of  surface-to-air  temperature  differences  within  an  enclosure.  Figure  8  represents 
the  experimentally  measured  deposition  velocity  to  the  vertical  surfaces  as  a  function  of  the 


Figure  8. 
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Surface-to-air  temperature  difference,  K 

Experimentally  measured  deposition  velocities  to  the  vertical  walls  of  a  1.2  x  1.2  x  1.2  m' 
enclosure  as  functions  of  particle  diameter  and  surface-to-air  temperature  difference 
(reprinted  with  permission  from  Thatcher  et  al 46  Deposition  velocities  predicted  theoreti¬ 
cally  by  Nazaroff  and  Cass65  for  a  vertical  isolated  flat  surface  are  also  plotted  for  compari¬ 
son.  Each  diamond  symbol  represents  the  average  value  obtained  by  10  local  surface 
extractions  over  the  central  portion  of  the  cool  or  the  warm  wall  from  a  single  experiment 
The  error  bars  span  the  maximum  and  minimum  values  of  measurements  at  each  location, 
while  an  arrow  is  used  to  indicate  that  the  lower  limit  of  the  95%  confidence  interval  for  the 
mean  was  less  than  the  detection  limit.  The  symbol  *  denotes  that  the  deposited  particle 
mass  was  under  the  detection  limit. 
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surface-to-air  temperature  difference  (±1.5K  and  ±10K.  respectively)  for  five  different  particle  sizes. 
The  deposition  velocity  predicted  theoretically  by  Nazaroff  and  Cass24h  was  also  plotted  for 
comparison. 

As  shown  in  Figure  8,  the  thermophoresis  effect  on  deposition  velocity  appears  to  be  diminishing  as 
particle  size  increases.  For  instance,  the  experimentally  obtained  deposition  velocity  to  the  cool  ver¬ 
tical  surface  (-10K)  for  0.1  pm  particles  is  more  than  3  orders  of  magnitude  higher  than  to  the  warm 
vertical  surface  (+10K).  On  the  other  hand,  the  measured  deposition  velocities  to  the  cool  and  warm 
walls  for  larger  particles  (1.3  and  2.5  pm)  are  approximately  the  same.  Furthermore,  for  particles  lar¬ 
ger  than  0.5  pm,  the  measured  deposition  velocity  on  the  warm  wall  at  +10K  is  higher  than  at  +  1.5K. 
This  observation  is  counterintuitive  since  a  warmer  surface  is  anticipated  to  enhance  thermophoretic 
repulsion  thus  reducing  particle  deposition. 

In  summary,  the  deposition  velocity  determined  experimentally  in  Thatcher  et  al.46  shows  relatively 
good  agreements  with  the  predictions  for  all  five  particle  sizes  studied  when  the  wall  surface  tem¬ 
perature  is  cooler  than  the  air.  In  contrast,  the  deposition  velocity  to  the  warm  wall  obtained  in  the 
experiments  appears  to  be  much  higher  than  predicted,  particularly  for  increasing  particle  size.  The 
discrepancy  between  experimental  findings  and  theoretical  predictions  may  be  attributed  to  the  fact 
that  the  existing  model  does  not  account  for  all  the  factors  influencing  particle  deposition  under  the 
settings  of  the  experiments.  The  surface-to-air  temperature  difference  influences  not  only  particle 
thermophoretic  velocity,  but  also  the  near-surface  airflow  pattern,  which  further  introduces  more 
complexities  toward  understanding  of  the  deposition  processes.  The  complicated  airflow  pattern  and 
flow  instability  in  their  chamber  experiments  are  believed  to  play  a  role  in  contributing  the  variations 
of  the  deposition  velocity  measurements.  Other  factors  such  as  particle  inertia  and  the  effect  of  cor¬ 
ners  on  airflow  not  being  addressed  in  the  model  may  be  of  importance  to  influence  the  deposition 
process. 


6.3.4  In  the  Presence  of  an  Electric  Field 

Electrostatic  forces  may  play  an  important  role  in  influencing  particle  deposition.  When  the  interior 
surface  of  an  enclosed  volume  has  the  tendency  to  acquire  electrostatic  charge  (e.g.,  made  of  poor 
conducting  materials),  a  local  near-surface  electric  field  may  develop,  which,  in  turn,  enhances  the 
deposition  of  charged  particles.*  Early  experiments  have  shown  that  the  presence  of  electric  charges 
on  surfaces  can  introduce  additional  variability  to  particle  deposition  measurements.66 

McMurry  and  Rader*  demonstrated  in  their  60-rn^  Teflon  chamber  experiments  under  natural  con¬ 
vection  that  the  enhanced  deposition  for  0.07-1  (im  particles  was  attributed  to  electrostatic  effects,  as 
seen  in  Figure  9.  In  another  set  of  Teflon  film  bag  (250  L)  experiments,  the  deposition  rates  of  singly 
charged  particles  could  be  significantly  increased  as  compared  to  particles  at  Boltzmann  distribu¬ 
tion,67  as  shown  by  the  symbols  V  and  T  in  Figure  10.  For  instance,  a  nearly  30  times  deposition 
enhancement  was  observed  for  the  0.1  (im  singly  charged  particles  in  comparison  to  particles  with 
Boltzmann  charge  distribution.  In  addition,  their  experimental  data  indicated  that  the  deposition  rates 


As  mentioned  in  Subsection  3.6,  airborne  particles  become  charged  by  collision  with  air  ions  and  carry  certain  charges 
depending  on  their  particle  size  according  to  Boltzmann  distribution 
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Deposition  rate  P,  s 


0.01  0.1  1.0  10.0 


Particle  diameter,  pm 

Figure  9.  Comparison  of  measured  particle  deposition  rates  in  a  60  m  Teflon  smog  chamber  with 
modeling  predictions  by  Crump  and  Seinfeld. AR  E  and  ke  denote  the  estimated  average 
electric  field  strength  on  the  chamber  interior  surface  and  the  parameter  of  airflow  tur¬ 
bulence  intensity,  respectively.  Note  that  the  experimental  data  exceeded  predicted  val¬ 
ues  for  0.07-1  pin  particles  when  electrostatic  effects  were  not  taken  into  account. 
Reprinted  with  permission  from  McMurry  and  Rader.67 


McMurry  and  Rader  [67]  Boltzmann  distribution  (n  =  0) 

McMurry  and  Radar  [67],  singly  charged  aerosols  (n  =  1) 

Shimada  at  al.  [69],  0  =  0 

Shimada  et  al  [69].  0>  =  50  V 

Shimada  at  al  [69],  0>  =  100  V 

Shimada  at  al.  [69],  0>  =  200  V 


Particle  diameter,  pm 

Figure  10.  Experimentally  determined  particle  deposition  rates  as  a  function  of  particle  diameter 
under  various  influences  of  electrostatic  attractions.  The  data  points  of  MuMurry  and 
Rader61  were  obtained  in  a  250-L  Teflon  film  bag  (natural  convection;  S/V  =  10.6  m"1),  in 
which  the  average  electric  field  was  estimated  as  45  V/cm  via  data  fitting  to  their  model. 
The  electrostatic  attractions  in  Shimada  et  al.69  were  established  by  applying  known  volt¬ 
age  (O)  at  the  bottom  of  the  enclosed  tank  (forced  convection;  S/V  =  47  nrf  r). 


31 


were  nearly  identical  for  oppositely  charged  particles  of  the  same  size.  The  laboratory  results  agreed 
satisfactorily  with  their  model  that  incorporated  electrostatic  drift  as  an  additional  transport  mecha¬ 
nism.  As  explained  in  the  model,  deposition  for  particles  larger  than  1  pm  and  smaller  than  0.05  pm 
were  still  dominated  by  gravitational  settling  and  Brownian/turbulent  diffusion,  respectively;  while 
electrostatic  effects  were  important  factors  influencing  deposition  of  0.05-1  pm  particles. 

Shimada  et  al.h<  performed  deposition  measurements  for  0.02-0.2  pm  particles  in  a  stirred  metal  tank 
where  the  turbulence  intensity  and  the  electric  field  could  be  controlled.  They  found  no  difference  in 
deposition  rates  for  charged  and  uncharged  particles  when  the  surfaces  were  grounded.  As  the  elec¬ 
tric  field  strength  increased,  however,  the  enhancements  of  particle  deposition  rates  were  observed,  as 
indicated  in  Figure  10. 

By  directly  collecting  particle  deposition  mass  onto  vertical  surfaces  in  a  1 .8-nr  cubic  chamber,  Lai  1 
reported  that  the  measured  particle  deposition  velocities  onto  acetate  sheets  for  3.5-9  pm  particles 
were  more  than  an  order  of  magnitude  higher  than  those  onto  glass  and  copper  plates,  owing  to  the 
electrostatic  effects.  The  laboratory  data  indicated  that  deposition  velocities  for  both  glass  and  copper 
surfaces  were  comparable,  and  the  use  of  anti-electrostatic  spray  on  the  surfaces  was  found  to  lead  to 
reduced  particle  depositions  by  minimizing  the  Coulombic  effect. 
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7.  Modeling  Particle  Deposition  and  Experimental  Validations 


Two  major  types  of  models,  both  theoretical-based  and  semi-empirical,  have  been  proposed  to  quan¬ 
tify,  predict,  and  explain  particle  deposition  in  enclosures.  The  first  type,  as  will  be  discussed  in  Sub¬ 
sections  7. 1  and  7.2,  involves  first  modeling  the  air  turbulent  structure  adjacent  to  the  enclosure  sur¬ 
faces,  and  particle  transport  is  subsequently  formulated  accounting  for  gravity,  diffusion,  and  other 
deposition  mechanisms.  One  key  challenge  in  this  modeling  approach  lies  in  the  determination  of  the 
near-surface  particle  eddy  diffusi vity ,  which  is  postulated  to  be  related  to  the  air  turbulence  intensity 
within  the  enclosure  as  well  as  the  distance  to  the  wall.  The  second  type  of  modeling  approach,  as 
will  be  described  in  Subsection  7.3,  applies  the  well-known  analogy  between  mass  and  heat  transfer, 
and  deposition  rate  of  particles  from  the  gas  phase  to  the  surface  due  to  diffusion  is  estimated  via  the 
mass  transfer  correlation. 

In  this  section,  the  modeling  developments  based  on  these  two  approaches  for  predicting  particle 
deposition  rates  in  enclosures  are  summarized,  together  with  the  available  experimental  investigations 
to  compare  against  the  modeling  analyses  as  well  as  the  limitations  associated  with  these  models. 


7.1  Homogeneous  Turbulence  Model 

Modeling  efforts  for  studying  particle  deposition  from  turbulent  flow  onto  enclosure  surfaces  were 
initiated  in  the  early  1950s.  Corner  and  Pendlebury^  developed  the  first  theoretical  model  accounting 
for  particle  deposition  onto  the  surfaces  of  horizontal  and  vertical  orientations  in  a  rectangular  enclo¬ 
sure  where  the  air  was  homogeneously  turbulent.  Their  model  was  derived  based  on  the  following 
key  assumptions:  (1)  the  air  outside  the  boundary  layer  is  homogeneously  turbulent,  and  aerosol  con¬ 
centration  is  uniform  (as  illustrated  in  Figure  1);  (2)  within  the  boundary  layer  of  thickness  S [  the 
fluid  motion  is  turbulent  with  random  fluctuations,  but  the  mean  fluid  motion  is  parallel  to  the  sur¬ 
face;  (3)  the  velocity  gradient  is  linear  within  the  boundary  layer;  (4)  the  mechanisms  of  gravitational 
settling,  Brownian/turbulent  diffusion  are  responsible  for  particle  transport  through  the  boundary 
layer;  (5)  turbulent  diffusion  dominates  Brownian  motion  at  the  edge  of  the  boundary  layer;  and  (6) 
particle  transport  is  quasi-steady-state  in  the  boundary  layer. 

Under  the  above  assumptions,  the  aerosol  concentration  in  the  boundary  layer  adjacent  to  the  surface 
is  governed  by 


dy 


-  vtsk  • 


ac 

dy 


=  0 


(21) 


where  k  is  the  unit  normal  vector  in  the  vertical  direction.  The  first  and  second  terms  of  Eq.  (21) 
account  for  the  processes  of  turbulent/Brownian  diffusion  and  gravitational  settling,  respectively. 
The  boundary  conditions  for  the  above  equation  are 
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(22) 


C  =  0  at  y  =  0 
C  =  Coo  at  y  >  8 

where  Coo  is  the  bulk  aerosol  concentration  outside  the  boundary  layer  The  particle  concentration 
profile  in  the  boundary  layer  can  thus  be  solved  from  Eqs.  (21)  and  (22).  In  the  Corner  and  Pendle- 
bury  model,  the  particle  eddy  diffusi  vity  is  approximated  using  Prandtl’s  mixing  length  theory  as 


Ep=key2, 


(23) 


where  is  the  turbulence  intensity  parameter  used  to  characterize  the  degree  of  turbulent  mixing 
inside  the  enclosure.  There  is  no  a  priori  way  to  estimate  the  values  of  k^  with  available  data;  hence, 
it  is  usually  determined  empirically  by  fitting  experimental  data  into  the  theory.  Comer  and  Pendle- 
bury53  suggested  that  k^  be  evaluated  as 


ke 


dU 
dy  ’ 


(24) 


where  K  is  the  von  Karmdn’s  constant  (usually  taken  as  0.4),  and  U  is  the  component  of  the  mean 
flow  velocity  parallel  to  the  surface.  The  velocity  gradient,  dU/dy,  is  approximated  by  means  of  fluid 
drag  force  balance  for  a  flat  plate  in  Corner  and  Pendlebury’s  model. 

Following  the  seminal  work  of  Comer  and  Pendlebury,^  several  additional  model  extensions  have 
been  proposed.  Crump  and  Seinfeld68  derived  expressions  for  estimating  particle  deposition  rate  in  an 
enclosure  of  arbitrary  shape,  and  the  analytical  solutions  for  a  spherical  vessel  were  specifically  for¬ 
mulated.  They  showed  that,  in  an  enclosure  having  only  horizontal  and  vertical  surfaces,  the  trans¬ 
port  processes  of  gravitational  settling  and  diffusion  can  be  treated  independently  by  vectorially 
summing  the  gravitational  settling  velocity  to  the  deposition  velocity  associated  with  the  diffusion 
process.  For  the  inclined  surfaces  of  a  sphere,  however,  their  derivations  indicated  that  these  two 
mechanisms  are  always  intimately  coupled  to  each  other  and  cannot  be  separated.  They  proposed  that 
the  value  of  k^  can  be  evaluated  from  the  fluid  energy  dissipation  rate,  instead  of  velocity  gradient  in 
the  boundary  layer  as  suggested  in  Comer  and  Pendlebury.^  Moreover,  Crump  and  Seinfeld68 
derived  the  particle  deposition  rate  with  a  more  general  form  of  £p  =  ke  yn,  where  n  can  be  any  num¬ 
ber,  which  could  be  obtained  by  empirically  fitting  experimental  data.  They  also  demonstrated  that 
the  use  of  the  exponent  n  =  3,  as  suggested  by  Friedlander,3  could  produce  analytical  expressions  to 
predict  deposition  that  are  analogous  to  those  with  n  =  2. 

In  summary,  Crump  and  Seinfeld^8  calculated  the  particle  deposition  velocity  onto  the  surface  of  an 
enclosure  to  be 
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where  v(s  is  the  particle  terminal  settling  velocity,  and  0  is  the  angle  between  normal  vector  of  wall 
and  gravity  direction  (rad). 

To  validate  their  own  theoretical  derivations,  Crump  et  al71  performed  particle  deposition  experi¬ 
ments  in  a  spherical  vessel.  By  fitting  the  data  and  using  n  =  2,  the  turbulence  intensity  parameter  kc 
in  their  study  was  evaluated  in  terms  of  the  volumetric  air  flowrate  into  the  chamber,  which  was 
related  to  the  turbulent  energy  dissipation  rate,  as  suggested  in  Okuyama  et  a  I.72  The  experimental 
results  in  Crump  et  al. 71  showed  good  agreement  with  the  analytical  expressions  with  respect  to  the 
dependence  of  particle  deposition  rates  on  particle  size  and  turbulence  intensity. 

Based  on  the  existing  model  developed  by  Crump  and  Seinfeld, 6s  McMurry  and  Rader"  incorporated 
electrostatic  effects  as  an  additional  particle  deposition  mechanism  to  evaluate  the  particle  deposition 
rate  in  an  enclosure.  Experiments  were  also  performed  to  measure  the  deposition  rate  of  neutral  and 
singly  charged  aerosols  in  a  250-L  Teflon  bag  and  a  60-m'  Teflon  smog  chamber,  respectively.  A 
good  fit  between  the  experimental  results  and  the  model  predictions  was  obtained  by  setting  n  =  2  and 
adjusting  the  values  for  ke  and  the  mean  electric  field 

Okuyama  et  al 23  studied  the  deposition  loss  of  monodisperse  aerosols  with  particle  diameters  of 
0.006  to  2  Jim  in  a  stirred  cylindrical  vessel.  The  turbulence  intensity  parameter  in  their  experi¬ 
ments  was  calculated  by  the  average  energy  dissipation  rate  per  unit  mass  of  air,  and  the  parameters 
of  k  and  n  were  obtained  by  fitting  the  experimental  data.  The  experimental  deposition  rates  com¬ 
pared  well  with  the  model  calculations  by  Crump  and  Seinfeld68  when  the  eddy  diffusivity  was 
assumed  to  be  proportional  to  the  2.7th  power  of  the  distance  from  the  surface  (i.e.,  n  =  2.7).  In  addi¬ 
tion,  more  deviations  from  the  model  were  observed  for  increasing  particle  size  or  flow  turbulence 
intensity,  and  this  is  likely  attributed  to  effects  of  enhanced  particle  inertia,  which  was  not  considered 
in  the  model. 

Cheng*'  measured  the  particle  deposition  rates  of  monodisperse  particles  ranging  from  0.005  to  2  pm 
in  a  spherical  chamber.  Both  the  chamber  temperature  and  air  velocity  profiles  under  various  turbu¬ 
lence  conditions  were  measured.  The  turbulence  intensity  in  the  experiments  evaluated  using  both 
methods,  including  the  velocity  gradient  as  shown  in  Eq.  (24)  and  energy  dissipation  rate,23, 72  gave 
reasonable  estimates  of  k<,  as  well  as  the  particle  deposition  rates.  The  data  were  well  explained  by 
Crump  and  Seinfeld’s  model,68  except  that  the  best-fitted  estimate  of  n  was  approximately  2.8,  instead 
of  n  =  2  as  supported  in  Crump  et  McMurry  and  Rader,67  and  Chen  et  al 42 


The  expression  ke  «  (elv'f12  suggested  in  Okuyama  et  al.12  originates  essentially  from  the  Prandd’s  mixing  length 

formula,  where  [el  v)^2  is  proportional  to  the  r.m.s.  velocity  gradient  in  the  boundary  layer.  The  symbol  £  refers  to  the 
turbulence  energy  dissipation  rate  and  v  is  the  kinematic  viscosity  of  the  fluid. 
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By  vectorially  adding  the  deposition  flux,  Nazaroff  and  Cass65  incorporated  the  thermophoresis  effect 
into  the  model  by  Corner  and  Pendlebury53  to  estimate  deposition  velocity  onto  a  vertical  isolated  flat 
plate  in  the  presence  of  surface-to-air  temperature  difference.  Their  analysis  showed  that,  as  particle 
size  increases,  an  increasingly  pronounced  difference  in  deposition  velocity  is  predicted  for  the  verti¬ 
cal  warm  relative  to  the  cool  walls.  The  experimental  findings  in  Thatcher  et  al .46  indicated  that  the 
model  appeared  to  greatly  underestimate  the  deposition  velocity  for  a  warm  surface,  especially  for 
supermicron  particles,  as  was  explored  in  Subsection  6.3.3. 

In  summary,  the  existing  experimental  data  agree  reasonably  well  with  the  theory  by  Corner  and  Pen- 
dlebury,53  including  its  extensions.  However,  there  are  still  some  issues  that  must  be  carefully  exam¬ 
ined  before  applying  the  model  to  various  scenarios.  For  example,  can  ke  and  n  be  determined  in  any 
particular  enclosure,  without  resorting  to  fitting  the  results  of  sophisticated  particle  deposition 
experiments?  These  two  parameters  are  key  and  vary  with  experimental  conditions;  they  ultimately 
influence  the  magnitude  of  the  particle  eddy  diffusivity  £p,  which  directly  affects  the  calculation  of 
deposition  rate  and  flux  to  the  surface. 

Table  4  summarizes  the  values  of  n  and  k^  determined  from  various  particle  deposition  experiments. 
The  scenario  of  n  =  2  is  regarded  as  the  “classical”  form  of  Corner  and  Pendlebury*  and  Crump  and 
Seinfeld.6*  This  is  also  supported  by  some  experimental  data  (Crump  et  al.;11  McMurry  and  Rader;67 
Chen  et  al.42).  In  other  experimental  findings,  however,  the  best  fit  of  the  data  to  the  model  occurs  at 
n  =  2. 6-2. 8  (Okuyama  et  al;  Van  Dingenen  et  al;  Holub  et  al ;  Cheng,  ),  close  to  n  =  3  as  sug¬ 
gested  by  Friedlander7  as  well  as  Pandian  and  Friedlander  4  based  on  a  theoretical  perspective  of  the 
analogy  between  mass  and  heat  transfer. 

>g  <  ••*21  ^  ..  _ | 

In  the  definition  of  particle  eddy  diffusivity  (in  cm  s  ,  £p  =  k<>y  )  when  n  =  2,  the  dimension  of  ke  (s  ) 

has  the  dimensions  of  a  rate  constant.  Non-integer  values  of  n  leads  to  Iq  with  a  dimension  of  L  "T  , 

which  not  only  lacks  a  solid  physical  foundation  but  also  causes  conceptual  and  practical  problems 

when  ke  is  to  be  evaluated  based  on  information  of  velocity  gradient  or  turbulent  energy  dissipation 

rate  (both  methods  assuming  n  =  2).  Based  on  the  rules  of  dimensional  analysis,  Benes  and  Holub 

suggested  a  modified  formulation  to  avoid  the  dimensional  inconsistency  problem: 


£p  =  ke5' 


(26) 


where  Sis  the  boundary  layer  thickness.  Using  this  new  expression  has  shown  to  yield  good  agree¬ 
ments  with  data  from  one  experimental  study. ^  However,  it  remains  unresolved  with  respect  to  the 
evaluation  method  for  n. 

Nevertheless,  Crump  and  Seinfeld6*  noted  that  the  choice  of  n  value  is  of  little  importance  from  the 
standpoint  of  their  theoretical  derivations.  Van  Dingenen  et  al.1  suggested  that,  after  re-examining 
the  data  by  Crump  et  al 11  and  McMurry  and  Rader, 6/  the  exact  value  of  n  is  trivial  as  long  as  an 
appropriate  k^  is  determined  in  an  independent  way,  for  instance,  the  evaluation  of  energy  dissipation 
rate  or  near-surface  velocity  gradient,  as  suggested  in  Okuyama  et  al 12 
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Table  4.  Summary  of  the  Values  of  n  and  ke  Determined  from  Various  Particle  Deposition  Experiments 

Particle  size 

Method  to  evaluate  K  Flow  mixing  method  Experimental  chamber  (lam)  Investigators 
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7.2  Three-Layer  Model 

Incorporating  the  information  on  the  structure  of  near-surface  turbulent  diffusivity,  Lai  and 
Nazaroff  5  proposed  a  model  in  which  only  one  parameter,  friction  velocity,  u\  is  required  to 
evaluate  particle  deposition  under  homogeneously  turbulent  conditions.  The  friction  velocity  is 
defined  by 


/  \ 

1/2 

( 

\ 

* 

dU 

u  = 

= 

v  — 

'  Pa  J 

dy 

O 

II 

>■* 

(27) 


where  Tw  is  the  shear  stress  at  the  surface,  pa  is  the  air  density,  and  v  is  the  kinematic  viscosity  of 
air.  The  parameter  u  is  intended  to  capture  the  turbulent  characteristics  in  the  vicinity  of  the  sur¬ 
face.  As  suggested  in  Eq.  (27),  the  friction  velocity  is  related  to  the  velocity  gradient  at  the  sur¬ 
face,  dU/dy,  which  can  be  evaluated  by  a  freestream  air  speed  IL  and  a  characteristic  length  of 
enclosure  surfaces  L:  6 


dU 

f  0.074^ 

PaVJ  1 

dy 

II 

O 

l  PaV  J 

2  J 

v  J 

(28) 


Another  approach  of  estimating  friction  velocity  is  to  measure  the  velocity  profile  in  the  loga¬ 
rithmic  flow  region  near  the  surface,  i.e.,  the  Clauser-plot  method.  Within  the  turbulent  bound¬ 
ary  layer,  the  time-averaged  velocity,  U,  as  a  function  of  distance  from  the  surface  y  is  expressed 
by 


2.5u*  J  yU^  > 
Uoo  v  j 


+  A , 


(29) 


where  A  is  a  constant.  Thus,  the  friction  velocity  can  be  inferred  from  the  slope  of  the  line  by 
plotting  measurements  of  U/IL  versus  the  logarithm  of  (yLL/v). 

Commonly  known  as  the  law  of  the  wall ,  a  turbulent  boundary  layer  consists  of  three  distinct 
zones  according  to  the  velocity  distribution  as  a  function  of  the  distance  perpendicular  to  the 
wall.76  The  approach  used  by  Lai  and  Nazaroff25  to  analyze  particle  deposition  from  turbulent 
flow  is  to  examine  the  turbulent  flow  structure  zone  by  zone,  and  to  formulate  particle  transport 
equations  for  each  zone,  assuming  (1)  gravitational  settling  and  Brownian  and  turbulent  diffusion 
are  responsible  for  the  particle  deposition  processes;  (2)  constant  particle  flux  in  the  concentration 
boundary  layer;  (3)  particle  eddy  diffusivity  well  represented  by  fluid  turbulent  viscosity;  and  (4) 
negligible  surface  roughness  effects.  Across  the  boundary  layer,  the  expression  for  the  deposition 
velocity  was  then  integrated  for  each  zone.  The  deposition  velocity  was  evaluated  for  vertical, 
upward,  and  downward  horizontal  surfaces,  and  the  first-order  deposition  rates  were  provided  for 
rectangular  and  spherical  cavities,  respectively.  This  approach  is  somewhat  complicated,  but 
remains  practical  to  use. 
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Table  5  summarizes  the  equations  required  for  calculating  particle  deposition  along  with  expla¬ 
nations  of  parameters  used  in  the  model.  Their  model  predictions  compared  well  with  published 
experimental  data  for  deposition  of  0.001-2  pm  particles  in  a  spherical  enclosure  (Cheng56),  as 
seen  in  Figure  1  1 .  Lai  and  Nazaroff25  indicated  that  their  model  calculations  yield  the  best  agree¬ 
ment  with  the  model  by  Crump  and  Seinfeld68  with  n  =  2.95,  very  close  to  n  =  3,  which  is  used  to 
characterize  convective  heat  transfer.74  It  is  expected  that  the  analogy  between  heat  and  mass 
transfer  should  hold  when  particle  inertia  is  insignificant. 


Table  5.  Summary  of  Equations  Used  for  Particle  Deposition  Analysis  in  the  Three-Layer  Model 
(reprinted  with  permission  from  Lai  and  Nazaroff25) 


Parameters 


Equations 


Integral  I 


I  =  j^3.64Sc^3(a  -  b)  +  39 J 


a  =— In 
2 


(,0. 


92Sc  1/3 +  4.3 


■r 


Sc"1 * *  +  0.0609 


-  tan"1 


8.6-I0.92SC"'/3 


•/t  I0.92SC  ’Z-’ 


1 

[l0.92Sc l/3  4  r+  j 

3 

+  >/3  tan  1 

2r+  -  10.92SC"1/3 

Sc 

+  7.669x  I0~4| 

3 

1  _ 

V?  I0.92SC"1/3 

Deposition  velocity 
verttcal  surface 


Deposition  velocity 

upward  horizontal  surface 


Deposition  velocity 

downward  horizontal  surface 


Deposition  rate, 

rectangular  enclosure 
Deposition  rate, 

spherical  enclosure1 


Nomenclature: 

Sc  -  v/D  =  parttcle  Schmidt  number 

v  =  kinematic  viscosity  of  air 

D  =  Brownian  diffusivity  of  particles 

t  =  dpu/2v 

dp  *=  particle  diameter 

u  =  friction  velocity 

vls  =  terminal  settling  velocity  of  particles 

Av~  area  of  vertical  surfaces 

Au  =  area  of  upward-facing  surfaces 

Ad=  area  of  downward-facing  surfaces 

V-  room  volume 

R  =  radius  of  spherical  enclosure 

D/A)  =  Debye  function  defined  by 


vd,u  - 


_ vjs 

f 

I -exp  - 


Vd,d  = 


Vts 


expl  -  ^  |  - 1 


h  vd  v^v  +  vd,u^u  +  vd,d^d 
P~  V 


Di 


dt 


v  The  integral  is  evelueted  anelyticalty  under  the  approximation  thet  Brownian  diffusivity,  D,  is  negligible  compared  with  eddy  diffusivity  for  y*  £  4.3,  where  y*  is  the 
normalized  distance  from  the  surface  This  epproximation  is  accurate  to  1%  or  better  for  particle  diameters  larger  than  0.01  pm  For  smaller  particles,  the 
integration  must  be  earned  out  numerically.  See  the  following  for  results 
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Particle  diameter,  dp  (pm) 

Integral,  1  (-) 

0.001 

29.1 

0.0015 

49.1 

0.002 

71.0 

0.003 

120.3 

0.004 

174.9 

0.005 

234.2 

0  006 

297.4 

0.007 

364.0 

0.008 

432.7 

0.009 

504.5 

0.01 

579.3 

*  In  the  limit  of  small  particles  (negligible  influence  of  gravitational  settling),  the  expression  simplified  to  p  =  3u*(RI)  In  the  limit  of  large  particles  (negligible 
influence  of  Brownian  diffusion),  the  expression  simplifies  to  p  =  3vtr(4R)  ’. 


Particle  diameter,  pm 

Figure  1 1 .  Predictions  of  particle  deposition  rate  as  a  function  of  particle  diameter  in  the 
three-layer  model  in  Lai  and  Nazaroff25  compared  with  experimental  data 
(Cheng56).  Reprinted  with  permission  from  Lai  and  Nazaroff.2^ 

Lai  and  Nazaroff48  performed  laboratory  experiments  to  measure  deposition  of  0.9-9  pm  particles 
from  turbulent  flow  onto  vertical  surfaces  of  a  cubic  aluminum  chamber.  The  experimentally 
measured  particle  deposition  velocities  for  9-pm  particles  were  higher  by  a  factor  of  30-1 50  as 
compared  to  their  own  model  predictions  for  friction  velocities  of  2.9-9. 8  m/s.  The  chamber 
walls  were  grounded;  thus,  electrostatic  force  was  considered  negligible.  They  note  that  the  dis¬ 
crepancy  between  experimental  observations  and  model  calculations  is  attributed  to  the  inadequa¬ 
cies  of  the  model,  in  which  the  key  transport  and  deposition  processes  for  supermicron  particles 
may  not  be  addressed  appropriately.  For  example,  particle  inertia  and  shear-induced  lift  force* 


A  particle  in  a  shear  flow  field  may  experience  a  lift  force  perpendicular  to  the  main  flow  direction  This  life  force 
arises  owing  to  particle  inertia  and  is  important  for  large  particles  Lai  and  Nazaroff45  postulated  that  shear-induced 
lift  force  may  be  important  in  their  experiments,  in  which  significant  velocity  gradient  adjacent  to  the  surface  is 
expected  owing  to  the  parallet  airflow  pattern  along  the  vertical  walls  of  the  chamber. 
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are  postulated  to  be  potentially  important  players  to  enhance  transport  of  larger  particles  through 
the  relatively  thin  boundary  layers  and  lead  to  subsequent  surface  deposition. 

Lai78  later  incorporated  the  mechanism  of  particle  inertia  into  the  three-layer  model  by  adopting 
the  commonly  accepted  electrical  resistance  analogy  in  atmospheric  dry  deposition  modeling,31  to 
simulate  supermicron  particle  deposition  from  turbulent  flow  to  a  vertical  surface.  Treating  the 
friction  velocity  as  the  fitting  parameter  and  incorporating  analyses  accounting  for  surface  rough¬ 
ness  effects,79  the  agreements  between  experimental  data  and  modeling  calculations  were  good 
for  smooth  surfaces,  but  less  satisfactory  for  rough  surfaces. 

The  influence  of  electrostatic  drift  was  considered  as  an  additional  particle  deposition  mechanism 
in  the  inertia-incorporated  three-layer  model  by  Chen  and  Lai.80  The  experimental  findings 
involved  with  electrostatic  effects  in  Lai70  again  demonstrated  that  the  particle  deposition  is  a 
complicated  process,  particularly  when  mechanisms  more  than  diffusion  and  gravitational  settling 
are  involved.  More  detailed  experiments  or  numerical  simulation  on  near-surface  airflow  struc¬ 
ture  will  be  helpful  to  shed  light  on  the  dynamics  of  particle  transport. 


7.3  Mass  Transfer  Model 

As  mentioned  in  Subsection  4.1,  particle  deposition  velocity  is  equivalent  to  a  mass  transfer  coef¬ 
ficient  with  the  units  of  length  per  time  (L'F1).  Based  on  the  analogies  for  mass  and  heat  trans¬ 
fer,81  the  correlations  for  transfer  coefficients  in  mass  and  heat  transport  should  be  of  the  same 
form  for  particle  mass  transfer  from  turbulent  flow  to  the  enclosure  surfaces  by  diffusion,  and  for 
convective  heat  transfer  in  jacketed  vessels,  assuming  that  the  average  roughness  heights  are 
immersed  within  the  viscous  sub-layer  (i.e.,  for  smooth  walls).  Following  this  analogy,  Pandian 
and  Friedlander74  proposed  a  semi-empirical  mass  transfer  expression  to  estimate  the  rate  of  parti¬ 
cle  deposition  due  to  turbulent  and  Brownian  diffusion  in  the  form  of  Sherwood-Reynolds- 
Schmidt  correlation: 


Sh=a(Re)2/3(Sc)1/'3, 


(30) 


where  Sh  is  the  Sherwood  number  =  kmDc/D,  km  is  the  mass  transfer  coefficient  or  particle  depo¬ 
sition  velocity  (m  s~l),  Dc  is  the  chamber  diameter  (m),  D  is  particle  diffusivity  (m2  s"1),  Re  is  the 
Reynolds  number  for  stirring  =  NDs/v,  N  is  the  impeller  speed  (rpm,  revolutions  per  minute)  in  a 
stirred  chamber,  Ds  is  the  stirrer  diameter  (m),  v  is  the  kinematic  viscosity  of  fluid  in  the  chamber 
(m2  s'1),  Sc  is  the  particle  Schmidt  number  =  v/D,  and  a  is  a  constant. 

By  fitting  the  experimental  data  in  Okuyama  et  al ,23  to  Eq.  (30),  Pandian  and  Friedlander74 
obtained  a  =  0.63  as  the  best  fit  for  particles  smaller  than  0. 1  pm  in  diameter  and  Re  no  greater 
than  3000.  Other  laboratory  results  were  also  used  to  compare  with  Eq.  (30),  and  good  agree¬ 
ments  were  found,  as  shown  in  Figure  12.  Using  different  sets  of  data,  the  values  of  fitted  a  have 
been  found  to  be  slightly  different. 

By  equating  the  expression  of  particle  deposition  rate  in  Crump  and  Seinfeld68  to  Eq.  (30), 
Cheng56  showed  that  the  mass  transfer  equation  proposed  by  Pandian  and  Friedlander  4  has  an 
equivalent  form  as  the  Crump  and  Seinfeld  model  with  n  =  3. 
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Figure  12.  Mass  transfer  correlation  of  particle  diffusional  deposition  in  smooth  walled 

stirred  chamber,  with  Sherwood  number  (Sh)  plotted  as  a  function  of  Re^Sc1^ 

The  dash  line  of  a  =  0.63  was  obtained  from  fitting  the  data  in  Okuyama  et 
a/.,23  and  the  straight  line  of  a  =  0.54  was  obtained  from  fitting  the  data  in  Oku¬ 
yama  et  alP  and  Cheng.56  Reprinted  with  permission  from  Cheng.56 

The  particle  deposition  rate  estimated  from  the  mass  transfer  correlation  thus  can  be  estimated  as 

|3  =  a  -5-  —  Re2/1  Sc1/3 , 

Dc  V 


where  S  and  V  are  the  interior  surface  area  and  the  volume  of  the  chamber,  respectively. 

Note  that  this  mass  transfer  correlation  indicated  in  Eq.  (30)  is  only  applicable  for  evaluating  par¬ 
ticle  deposition  rate  in  the  diffusional  regime  (e.g.,  dp  <  0.1  Jim)  in  small  cylindrical  or  spherical 
vessels  with  smooth  surfaces,  assuming  that  Re  can  be  defined  as  stated  in  Eq.  (30).  For  larger 
particles  and  higher  fluid  turbulence,  deviations  from  Eq.  (30)  were  observed  because  diffusion 
was  no  longer  the  only  mechanism  for  deposition  and  other  processes,  such  as  inertial  drift,  may 
contribute  to  particle  transport  to  the  surfaces. 


(31) 
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8.  Summary 


The  phenomenon  of  particle  deposition  onto  surfaces  in  a  confined  environment  is  frequently 
encountered  in  numerous  industrial  and  environmental  settings,  either  as  a  desirable  or  an  unde¬ 
sirable  outcome.  This  report  presented  a  literature  review  on  the  important  physical  factors  that 
influence  particle  deposition  in  enclosures,  as  well  as  the  available  experimental  techniques  and 
modeling  approaches  used  for  characterizing  the  rate  of  particle  deposition.  Experimental  find¬ 
ings  compared  with  the  model  calculations  were  presented,  and  the  caveats  with  respect  to  the 
models  were  discussed. 

As  was  already  known,  and  is  further  substantiated  by  the  scientific  evidence  presented  here,  the 
process  of  particle  deposition  is  a  complicated  phenomenon.  Transport  behavior  of  particles  to 
surfaces  is  governed  by  the  nature  of  airflow  near  surfaces,  including  the  turbulence  intensity,  the 
surface  characteristics,  and,  importantly,  particle  size.  The  experimental  studies  reviewed  in  this 
report  have  revealed  that  the  deposition  rate  varies  broadly  across  conditions,  and  the  measure¬ 
ment  results  are  affected  by  various  factors  acting  simultaneously.  Direct  measurements  of  parti¬ 
cle  deposition  onto  surfaces  of  interest  are  a  challenging  task  to  perform,  and  the  contribution 
from  various  potential  deposition  mechanisms  under  realistic  circumstances  makes  modeling 
work  difficult.  Nevertheless,  further  progress  to  elucidate  the  processes  of  particle  deposition 
will  require  continued  efforts  to  conduct  more  carefully  controlled  experimental  investigations,  in 
which  the  particle  size,  near-surface  air  flow  conditions,  and  the  nature  of  surfaces  are  well 
characterized. 
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PHYSICAL  SCIENCES  LABORATORIES 


The  Aerospace  Corporation  functions  as  an  "architect-engineer"  for  national  security  programs,  specializing  in 
advanced  military  space  systems.  The  Corporation’s  Physical  Sciences  Laboratories  support  the  effective  and 
timely  development  and  operation  of  national  security  systems  through  scientific  research  and  the  application  of 
advanced  technology.  Vital  to  the  success  of  the  Corporation  is  the  technical  staff  s  wide-ranging  expertise  and 
its  ability  to  stay  abreast  of  new  technological  developments  and  program  support  issues  associated  with  rapidly 
evolving  space  systems.  Contributing  capabilities  are  provided  by  these  individual  organizations: 

Electronics  and  Photonics  Laboratory:  Microelectronics,  VLSI  reliability,  failure  analysis, 
solid-state  device  physics,  compound  semiconductors,  radiation  effects,  infrared  and  CCD 
detector  devices,  data  storage  and  display  technologies;  lasers  and  electro-optics,  solid-state 
laser  design,  micro-optics,  optical  communications,  and  fiber-optic  sensors;  atomic  frequency 
standards,  applied  laser  spectroscopy,  laser  chemistry,  atmospheric  propagation  and  beam 
control,  LIDAR/LADAR  remote  sensing;  solar  cell  and  array  testing  and  evaluation,  battery 
electrochemistry,  battery  testing  and  evaluation. 

Space  Materials  Laboratory:  Evaluation  and  characterizations  of  new  materials  and 
processing  techniques:  metals,  alloys,  ceramics,  polymers,  thin  films,  and  composites; 
development  of  advanced  deposition  processes;  nondestructive  evaluation,  component  failure 
analysis  and  reliability;  structural  mechanics,  fracture  mechanics,  and  stress  corrosion;  analysis 
and  evaluation  of  materials  at  cryogenic  and  elevated  temperatures;  launch  vehicle  fluid 
mechanics,  heat  transfer  and  fight  dynamics;  aerothermodynamics;  chemical  and  electric 
propulsion;  environmental  chemistry;  combustion  processes;  space  environment  effects  on 
materials,  hardening  and  vulnerability  assessment;  contamination,  thermal  and  structural 
control;  lubrication  and  surface  phenomena.  Microelectromechanical  systems  (MEMS)  for 
space  applications;  laser  micromachining;  laser-surface  physical  and  chemical  interactions; 
micropropulsion,  micro-  and  nanosatellite  mission  analysis;  intelligent  microinstruments  for 
monitoring  space  and  launch  system  environments. 

Space  Science  Applications  Laboratory:  Magnetosphenc,  auroral  and  cosmic-ray  physics, 
wave-particle  interactions,  magnetospheric  plasma  waves;  atmospheric  and  ionospheric  physics, 
density  and  composition  of  the  upper  atmosphere,  remote  sensing  using  atmospheric  radiation; 
solar  physics,  infrared  astronomy,  infrared  signature  analysis;  infrared  surveillance,  imaging  and 
remote  sensing;  multispectral  and  hyperspectral  sensor  development;  data  analysis  and 
algorithm  development;  applications  of  multispectral  and  hyperspectral  imagery  to  defense,  civil 
space,  commercial,  and  environmental  missions;  effects  of  solar  activity,  magnetic  storms  and 
nuclear  explosions  on  the  Earth’s  atmosphere,  ionosphere  and  magnetosphere;  effects  of 
electromagnetic  and  particulate  radiations  on  space  systems;  space  instrumentation,  design, 
fabrication  and  test;  environmental  chemistry,  trace  detection;  atmospheric  chemical  reactions, 
atmospheric  optics,  light  scattering,  state -specific  chemical  reactions,  and  radiative  signatures  of 
missile  plumes. 
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